Berry phases and pairing symmetry in Holstein-Hubbard polaron systems 
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We study the tunneling dynamics of dopant-induced hole 
polarons which are self-localized by electron-phonon cou- 
pling in a two-dimensional antiferromagnet. Our treatment 
is based on a path integral formulation of the adiabatic 
(Born-Oppenheimer) approximation, combined with many- 
body tight-binding, instanton, constrained lattice dynamics, 
and many-body exact diagonalization techniques. The ap- 
plicability and limitations of the adiabatic approximation in 
polaron tunneling problems are discussed in detail and adi- 
abatic results are compared to exact numerical results for a 
two-site polaron problem. Our results are mainly based on 
the Holstein-t J and, for comparison, on the Holstein-Hubbard 
model. We also study the effects of 2nd neighbor hopping and 
long-range electron-electron Coulomb repulsion. The polaron 
tunneling dynamics is mapped onto an effective low-energy 
Hamiltonian which takes the form of a fermion tight-binding 
model with occupancy dependent, predominantly 2nd and 3rd 
neighbor tunneling matrix elements, excluded double occu- 
pancy, and an effective intersite charge interactions. Anti- 
ferromagnetic spin correlations in the original many-electron 
Hamiltonian are reflected by an attractive contribution to the 
1st neighbor charge interaction and by Berry phase factors 
which determine the signs of effective polaron tunneling ma- 
trix elements. In the two-polaron case, these phase factors 
lead to polaron pair wave functions of either d^2_y2-wa.Ye sym- 
metry or p-wave symmetry with zero and nonzero total pair 
momentum, respectively. Implications for the doping depen- 
dent isotope effect, pseudo-gap and Tc of a superconducting 
polaron pair condensate are discussed and compared to ob- 
served properties of the cuprate high-Tc materials. 

74.72.-h, 71.38.-Hi, 75.10.Lp, 71.27.-Ha 



I. INTRODUCTION 

The symmetry of the superconducting order pa- 
rameter in the cuprate high- Tc superconductors had 
been controversial 0| before phase-sensitive experiments 
firmly established the d^2_y2-waYe pairing symmetry 
in YBa2Cu307, using tricrystal ring magnetometry 
SQUID interferometry and single-junction modula- 
tion Q. Migdal-Eliashberg-type diagrammatic theories 
find d-wave pairing to be favored by antiferromagnetic 
(AF) spin-fluctuation exchange ||^ and s-wave pairing by 
the conventional electron-phonon mechanism, jsj]^ There 
is indeed strong experimental evidence for the impor- 



tance of both AF spm 
phonon interactions 



correlations and electron- 
in the cuprates. However, when 
combined in the diagrammatic approach, the two mech- 
anisms are mutually destructive, since d-wave pairing is 
strongly suppressed by phonons and s-wave pairing is 
suppressed by AF spin fluctuations, respectively. Also, 
the magnitude of the observed isotope effect in cuprate 
systems away from "optimal" doping points towards 
an unusually strong electron-phonon effect which cannot 
be accounted for in the diagrammatic approaches. |^ 

Strong-coupling studies, |lO|"[l^ going beyond the 
Migdal-Eliashberg regime, suggest that the AF spin cor- 
relations themselves can effectively enhance the electron- 
phonon effect, by lowering the electron-phonon coupling 
threshold for polaron formation, that is, the threshold 
for electron-phonon induced self-localization ||l^ of the 
dopant-induced carriers in the Cu02 planes. In the 
present paper, we show how the tunneling dynamics of 
such self-localized holes in an AF correlated spin back- 
ground may lead to d- and other non-s-wave pairing 
states which are not suppressed by coupling to the lattice 
degrees of freedom. 

A Berry phase factor in finite systems with time- 
reversal symmetry has been relevant to the observation of 
half-odd-integer quantum numbers in the spectrum of the 
Na3 molecule, , to the cross section of the H-I-H2 re- 
action and its isotope analogs, ||l^ and to the problem of 
integer vs. half-odd-integer spin tunneling in anisotropic 
potentials. [0 Contributions to the pair binding energy 
in the Cgo molecule have also been discussed in terms 
of Berry phase arguments. [|8| In the present case, the 
non-s-wave symmetry is caused by a (—1) Berry phase 
factor, associated with predominantly second- and third- 
neighbor polaron tunneling processes. It also determines 
the total momentum: the one-polaron ground state has a 
momentum on the Fermi surface of the half-filled tight- 
binding model on the square lattice. The dynamics of 
few hole polarons reflects the local AF spin correlations 
of many electrons through the Berry phase factor. 

This paper is organized as follows: In Sec. ||, we intro- 
duce the basic Holstein-Hubbard and Holstein-t J model 
Hamiltonians, and their extensions to include 2nd neigh- 
bor hopping or long-range Coulomb repulsion. We then 
derive the effective action for the lattice degrees of free- 
dom in the adiabatic approximation. In Sec. Ill, we illus- 



trate the basic physical principles and formal concepts of 
our adiabatic treatment of the polaron tunneling in the 



1 



context of a simple two-site model. In Sec. IV, we discuss 



the conditions under which the adiabatic approximation 
is valid, as well as its limitations when applied to pola- 
ronic systems on large / macroscopic lattice systems. In 
particular, we clear up some recent misunderstandings 
concerning the applicability of the adiabatic approach to 
polaronic systems. In Sec. we use an instanton ap- 
proach to elucidate the basic structure of the low-energy 
tunneling dynamics of hole polarons in Holstein-Hubbard 
or Holstein-t J systems near half-filling. We show that the 
dynamics of such hole polarons is governed by an effective 
tight-binding Hamiltonian which includes 2nd and 3rd 
neighbor hopping matrix elements and a 1st neighbor at- 
traction. In Sec. VI, we discuss the Berry phase factors 
and, with the help of lattice symmetry operations, we 
show how such phases can be properly assigned to each 
segment of a closed tunneling path. The Berry phase fac- 
tors are then interpreted in terms a quasiparticle statis- 
tics and internal symmetries of the many-electron wave 
functions. In Sec. VII, we analytically solve the effective 



model to show how the Berry phase factors determine the 
total momenta and internal symmetries of the few-hole- 
polaron wave functions. In Sec. VIII, we report numeri- 



cal results for the effective polaron hopping and effective 
pair binding energy as functions of the phonon frequency 



and electron-phonon coupling strength. In Sec. [X 



discuss the implications of our numerical results for a 
possible superconducting pairing instability, the isotope 
effect and the pseudo-gap in a hole polaron liquid at finite 
doping concentration in the nearly half-filled Holstein-i J 
and -Hubbard systems and compare to experimental ob- 
servations in the cuprates. In Sec. we summarize the 
present work. Part of the results presented in this pa- 
per were reported briefly in an unpublished paper and 
proceedings. [H 



II. MODEL AND EFFECTIVE ACTION 



We use mainly the Holstein-t J model [^,^ and oc- 
casionally the Holstein-Hubbard model for comparison. 
Later, we also include 2nd neighbor electron hopping 
and/or long-range electron-electron repulsion terms in 
the model. The total Hamiltonian is of the general form 



H — He + He 



-ph 



H 



ph 



(1) 



where H^ is the purely electronic i J or Hubbard model 
part, defined on a two-dimensional (2D) square lattice 
with lattice sites j = 1 . . .N and on-site electron occupa- 
tion numbers specified below. 



He- 



ph 



(2) 



is the Holstein electron-phonon (EP) interaction, cou- 
pling the local oscillator displacement Uj to the electron 
on-site occupation rij with an EP coupling constant C 
and 



H 



M 



(3) 



describes the non-interacting Einstein phonon system, 
consisting of the bare harmonic lattice potential H^, with 
restoring force constant K, and of the lattice kinetic en- 
ergy Hm with an atomic mass M and conjugate momenta 
Pj = —ihd/duj. If we rescale to dimensionless displace- 
ments and conjugate momenta 



Uj = Uj/u-p, Pj = —id/duj , 
with the small polaron shift 

C 

Up = — , 
K ' 



(4) 



(5) 



then -ffo-ph and i?ph can be completely parametrized in 
terms of only two characteristic energies, the bare Ein- 
stein phonon energy 



(6) 



and the ionic-limit [t — > 0) small polaron binding energy 



Ep = — . 
K 



(7) 



All results in the following are therefore stated in terms 
of wp, 17 and Ep only ||ll|,|l|l. 
The tJ model is written as l20f| 



He = - 



H.c, 



J 



riiUj • 



(8) 



with 1st neighbor electron hopping t and AF exchange 
coupling J. Here, Cirr annihilates an electron with 



spin a at site i, rii^ — c[^c. 



\ J2a,i3 cla'^apCip with (7 = {ax,(Jy,(Jz) denoting the vec- 
tor of Pauli spin matrices. The Hilbert space is restricted 
to states with no double occupancy at any site j, i.e., 
rij = 0, 1 only. 

The Hubbard model is written as 



He = -t 



{id),"- 



H.c] +C/^n,Tn,x (9) 



with on-site repulsion U and no restrictions on the on- 
site occupancy, i.e., rij = 0, 1, 2. In the following, we set 
h = 1, t = 1 and use J = 0.5t or U = 8t in the tJ or 
Hubbard model, respectively, unless stated otherwise. 

In addition to the standard tJ and Hubbard electronic 
model, we will also study the effects of additional, po- 
tentially important electronic terms, the 2nd neighbor 
hopping Hf and the long-range Coulomb repulsion Hic- 
Namely, 
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(10) 



where {i,j} denotes 2nd neighbor bonds and t' is the 
corresponding 2nd neighbor matrix element. The long- 
range 1/r Coulomb repulsion is 



(11) 



where i and j are summed independently over all sites 
excluding i = j and denotes the vector pointing from 
i to j, measured in units of the 2D lattice constant a = 
1. On a lattice with periodic boundary conditions we 
make the definition of \rij \ unique by requiring to be 
a vector of the shortest possible length connecting i to 
j, subject to all possible periodic boundary shifts. The 
matrix element Vc is thus the Coulomb repulsion energy 
between two electrons at 1st neighbor distance. 

To study the tunneling dynamics of self- localized holes, 
we consider the path integrals for transition amplitudes 
in imaginary time in the Born-Oppenheimer (adiabatic) 
approximation. Following the standard Feynman- Trotter 
approach pl| , we break up the Hamiltonian in the 
imaginary-time evolution operator 



lim I e 

— *oo \ 



-ArHo 



(12) 



where At = (3/L, Hm is the lattice kinetic energy defined 
in Eq. (^) and the 0th order part Ho = H — Hm com- 
mutes with all lattice displacement operators Uj . At each 
time slice = fcAr, with k = 1 . . . L, we now insert a 
complete set of electron-phonon basis states | xl''') which 
are chosen to be simultaneous eigenstates of H^ and of 
all Uj . They can be written in the form 



(13) 



where | is the lattice part and | ^''■''^(u)) the elec- 
tronic part of I Xm*^'')- Written in 1st quantized notation, 
the lattice part is simply 



^u{x) = 5{u-x) 



(14) 



with lattice coordinate vectors x = {xi . . .xn) and u = 
{ui . . . Mat). The electronic part | ^''^'''(w)) denotes the k- 
th electronic eigenstate of the 0th order adiabatic Hamil- 
tonian 



Ho{u) = H^ + i?o-ph(w) + Hk{u) , 



(15) 



at fixed u. That is, Ho{u) is defined to act only on the 
electronic degrees of freedom at fixed (c-number) lattice 
displacement coordinates u = (ui . . . un) and 

Hoiu) I *(")(«)) = W^''\u) I *('^'(u)) (16) 



where | ^'(")('')) and its eigenenergy Wq'^\u) depend 
parametrically on the lattice displacements u. The ex- 
act imaginary time evolution under H can thus be repre- 
sented by a path integral with a Euclidean action, written 
at finite L as 



L 

5Kr),«(r)]=5][(Af/2) 



fc=i 



At 



+ ATWt''\uiTk)) 

-ln(vl/(«'=)(«(r,)) I M/(«'=-^)(«(t,_i)))] ■ (17) 

The path integration is to be carried out both 
over the continuous lattice coordinates u(t/j) = 
(Mi(rfc) . . . UN{Tk)) and over the discrete electronic quan- 
tum numbers Kk = n{Tk). 

In the 0th order adiabatic approximation, correspond- 
ing formally to the M —t oo limit, one neglects the imag- 
inary time evolution of u altogether and replaces uu by a 
T-independent classical field. The 1st order adiabatic ap- 
proximation restores the r-dependence of the lattice co- 
ordinates u, under the simplifying assumption that the 
electrons follow the motion of the lattice adiabatically. 
That is, the path integration is restricted to configura- 
tions where, during r-evolution, the electrons remain in 
the same eigenstate, i.e., k/j = Kfc-i = k = const. Transi- 
tions between different electronic eigenstates ^ Hk-i 
are neglected. Formally, this approximation restores the 
leading order 1/M corrections to the lattice dynamics. 
At sufficiently low temperatures, one restricts the path 
integral further to include only the electronic ground- 
state K — 0. Suppressing the (K)-superscript altogether, 
one then arrives at the standard 1st order adiabatic 
(Born-Oppenheimer) approximation, with an effective 
Euclidean action 



^A.Kr)]=E[fE^ 



,-(Tfe) -Uj(Tfc_l)]2 



At 

+AtWo{u{tu)) 
-H^{u{Tk)) I ^'(u(rfc_i) 



(18) 



Note that Sxd depends explicitly only on the it- 
coordinates of the lattice. The first (M/2-) term is the 
standard form of the lattice kinetic energy for discretized 
imaginary time (finite L). The electronic groundstate en- 
ergy Wo(u) plays the role of a Oth-order (in 1/M) effec- 
tive lattice potential energy. The last term, containing 
the logarithms of the electronic groundstate wavefunc- 
tion overlaps at adjacent time slices during r-evolution, 
contains the Berry phase and 1 /M corrections to the po- 
tential energy, as we will now discuss. 
In exp(— S'Ad[i*(T)]), the overlap product 

L 

Q[u{t)] ^ W (f (u(rfe)) I ^>(u(Tk-i))) (19) 
fe=i 

enters which contains the Berry phase factor. 
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exp(-z0[«(T)]) = Q[u{t)]I I Q[u{t)] 



(20) 



i.e., 6'[m(t)] — — Imln(Q[it(T)]). Due to time-reversal 
symmetry, all | ^'(u(rfc))) have real amplitudes in an ap- 
propriately chosen electron basis and hence the phase fac- 
tor is real: exp(— ^^[-^(t)]) = sig\i{Q[u{T)\) . Taking L 
oo, we can also rewrite Re(ln (5[m(t)]) = In | Q[u(r)] | 
in S'AdlMC'r)] as a 1/M correction to the effective lattice 
potential which thus becomes 



W{u) = Wo{u) + Wi(u) , 
with Wi given by 



Thus, the effective action for L ^ oo becomes 



(21) 



(22) 



k=l ] 

+/^tW{u{tu))\ +ie[u{T)] 



(23) 



Equivalent results can be derived in the Hamiltonian 
approach to the adiabatic approximation. The basic idea 
here is to restrict the full electron-lattice Hilbert space 
to an "adiabatic" subspace which is spanned by the set 
of 0th order adiabatic electron-lattice eigenstates | Xu^) 
defined above in Eqs. (p^)-(p^) with k restricted to the 
electronic groundstate n — Q. The adiabatic subspace 
thus consists of EP states of the general form 



(24) 



where 4'i'u) is an arbitrary (square-integrable) wavefunc- 
tion amplitude which depends only on the lattice coordi- 
nates u. The basic approximation step is then to project 
the full EP Hamiltonian H onto the adiabatic subspace. 
In this manner one arrives at a 1st order effective Hamil- 
tonian H^d which is mathematically equivalent to the 
1st order adiabatic Euclidean action S'Ad in (p3|), after 
L — > oo. Since the adiabatic EP states | (f>) can be ex- 
pressed entirely in terms of their wavefunction amplitude 
4>{u), one can recast if Ad into the form of an effective 
Hamiltonian acting only on the lattice coordinates u in 
4>{u), without explicit reference to the underlying elec- 
tronic groundstate wavefunction | ^'1*''' ) contained in | 0) . 
However, it is crucial to keep in mind the formal relation- 
ship (^4|) between the full adiabatic EP state | (j)) and its 
lattice wavefunction amplitude 4'{u) if one wants to prop- 
erly compare 1st order adiabatic results to exact results, 
obtained by e.g. numerically diagonalizing the full EP 
Hamiltonian on small model clusters . 

In systems obeying standard harmonic lattice dynam- 
ics, the 0th order Born-Oppcnheimer "energy surface" 
Wo(u) exhibits one unique global minimum configura- 
tion which is, in terms of energy or in terms of 



configurational distance, well separated from other, if ex- 
istent, local minima. In that case, the path integral is 
dominated by small-amplitude "harmonic" fluctuations 
around and a description of the lattice dynam- 

ics in terms of renormalized harmonic oscillators, i.e., 
phonons, remains valid. Since the displacement excur- 
sions around 

y(min) 

small, so are the fluctuations 
in the electronic wavefunction | ^'(u)); hence the small- 
amplitude ("phonon") paths all have 6'[M(r)] — and 
Berry phase effects are negligible. Also, the u-derivatives 
of I ^'(u)) entering into Wi are well behaved and the 
TO-th order u-derivatives of the overlap matrix elements 
{duj^{u) I duj'i'iu)) are typical of the order of inverse 
lattice constants or inverse atomic distances raised to the 
(to -|- 2)-th power. The M^i -contribution to the harmonic 
restoring force constants, for example, are thus smaller 
than the 0th order Wo-contributions by factors of order 
of the 4th power of the lattice oscillator zero-point dis- 
placement amplitude over the lattice constant. Thus, the 
electronic overlap factor effects Wi{u) and 6'[u(r)] can be 
altogether neglected. 

By contrast, in polaronic systems the 0th order lattice 
potential Wq exhibits a large number of nearly degenerate 
local minima. The low-energy lattice dynamics is dom- 
inated by tunneling processes between the local minima 
which requires anharmonic large-amplitude excursions of 
the local displacement coordinates Uj and large local re- 
arrangements of the electronic wavefunction | ^'(u)) [0. 
In that case electronic overlap effects arising from both 
6'[u(t)] and Wi{u) can become quite important. 



III. TWO-SITE PROBLEM 

The two-site version of the Holstein-Hubbard model 
|p3| , p4| (Q-jl) is a simple toy problem which retains some 
essential physical features of the lattice polaron prob- 
lem. We will use it here to elucidate the basic underlying 
physical ideas and formal concepts of our adiabatic treat- 
ment of polaron formation and polaron tunneling dynam- 
ics and, also, to test the validity and illustrate some im- 
portant limitations of the adiabatic approximation. We 
restrict ourselves to the single-electron case on two sites, 
with an electronic intersite hybridization t. Hence there 
are no correlation (Hubbard-C/) effects and the adiabatic 
electronic wavefunction | vE'(m)) can be solved exactly by 
diagonalizing Hq[u) which reduces to a 2 x 2 matrix. 

The two sites are labeled 1 and 2 with on-site oscilla- 
tor coordinates ui and U2 and on-site electron occupation 
numbers ni and 712. Introducing symmetrized displace- 
ment coordinates 



u± = {ui ± U2)/\/2 , 

Wq and Wi can be written as 



Wo(u) = iyo+(u+) + W^o-(u-) , 



(25) 



(26) 



where 
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K C 

Wo+{u+) = —u\ + 

2 



1 

2 Vup 



\/2 Vup 



(27) 



2V^y Vsv^ 



iJp, (28) 



and 



Wiiu) = Wi-iu^) = - 



1 



4:Ep 



(29) 



There is no Berry phase, i.e., 9[u{t)] = Q, and the problem 
becomes equivalent to solving the Hamiltonian dynamics 
of a (fictitious) quantum particle of mass M moving in 
a two-dimensional (u+, u_)-plane subject to the effective 
potential W{u) = Wo{u) + Wi{u). Because of (26) and 
( p9| ) this dynamics is separable when written in terms of 
u+- and u_-coordinates. 

Since u+ couples only to the total electron charge n+ = 
ni + n2, which is conserved, the M^o+-part of Wq is just a 
harmonic potential, with its equilibrium position shifted 
to 



,(0) 



-Up 



/V2 



(30) 



by the constant pulling force Crt+ /\/2 exerted by the to- 
tal electron charge. Also, the electron groundstate wave- 
function I ^'(w)) =1 5'(m_)), and hence Wi, does not 
depend on u+. The dynamics of the u+-coordinate is 
therefore trivial, at least for processes conserving the to- 
tal electron number. 

Since m_ couples to the charge imbalance rt_ = ?ii — 712 
between the two sites, the shape of Wq- is renormal- 
ized by the EP coupling and Wi contributes to the u_- 
dynamics. We first consider the 0th order contribution 
Wo- , shown in Fig. |^(a) for several iJp-values. 

For small Ep, Wq- retains a single global minimum at 
U- = 0. Its harmonic restoring force 



■H^o-(.-^0)^(l-f)ir 



(31) 



softens with increasing Ep and changes sign when Ep 
reaches a critical value 



E, 



(crit) 



= 2t 



(32) 



For Ep > Ep''^*'\ the character of Wq- changes quali- 
tatively: Wo- acquires two degenerate minima at U- — 
±u^', separated by a maximum at u_ = 0, with 



,(0) 



/ 2t \2 up 



(33) 



where up = C/K is the polaron shift (^. w'^'' approaches 
up/\/2 in the strong-coupling limit 



Ep > t 



(34) 



The height of the 0th order potential barrier separating 
the two minima, 



A 



so 



Wo-iO) - Wo-{±u^^^) - (i - -^Yep (35) 



increases with Ep and approaches Ep /4 in the strong- 
coupling limit (34). 



The physical origin of the double-well potential can be 
most easily understood starting from the "ionic" {t = 0) 
limit of the model: For t — 0, the two electronic eigen- 

states of Hq{u), 



(«)\ = 



and 



= 1 ni = 1, n2 = 0) 



ni 



0, 712 = 1) 



(36) 



(37) 



have the electron completely localized on site 1 and 2, 
respectively, with eigenenergies 



± up/V2f - ^Ep , (38) 



where the upper (lower) sign refers to {W'^'^^K as 

shown by the two parabolic potential curves in Fig. 0(a). 
Assuming C > 0, | vl/(^i)) is the groundstate for u- < 
and I vl/(^2)^ ^ Q y_ _ ^Yie two parabolic 

eigenenergy curves W^(")(u_) and W'^^'^\u-) intersect, 
both states are degenerate and the groundstate wavefunc- 
tion changes discontinuously as a function of u_. When 
the hybridization t is turned on, the two fully localized 
wavefunctions | 'i'^^^)) and | v[/(^2)^ become mixed, the 
electronic degeneracy at u_ = is lifted and a mini- 
mum excitation gap of 2t opens up between the electronic 
ground- and 1st excited states. The sharp cusp at m_ = 
in the t = double-parabolic potential function 



Wo-iu-) |t=o -min(T^(")(u_),W'(^2)(^_)) 



yd- 



-Up 



(39) 



[see Fig. |l|(a)] is rounded by the finite t; as a function of 
M_, the groundstate wavefunction | ^'(u_)) now changes 
continuously at u_ = 0. However, | ^'(u_)) still has 
predominantly | ^'(^^^)-character near u_ ~ — u^-* and 

predominantly | '5 ^^2^) -character near m_ — u^-*. With 
increasing t, the tunneling barrier height decreases, ini- 
tially by about t. The barrier vanishes when t reaches 
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^(crit) _ which is equivalent to the above condition 

(H), for 4"'*\ 

If one examines the groundstate wavefunction | 
and its charge distribution | rij \ "^{uJ)) for 

near the two potential minima ±u^^\ one 
thus finds the electron predominantly localized at site 
1 when 7i_ = — w'^'' and predominantly at site 2 when 

M_ = +u'i^\ assuming again C > here and in the fol- 
lowing. By contrast, at the potential minimum — 
in the regime Ep < Ep"*\ the electron charge is delo- 
calized evenly between the two sites 1 and 2. Thus, at 
the level of the 0th order adiabatic approximation (i.e., 
neglecting the lattice kinetic energy), the transition from 

the single-well potential case Ep < Ep to the double- 
well case Ep > Ep"^"^ is essentially a transition from a 
delocalized non-degenerate groundstate (m_ = 0) to a lo- 
calized degenerate groundstate (m_ — ±u^''). In the for- 
mer case, the electron's delocalization energy dominates 
and it is energetically favorable for the electron wavefunc- 
tion to be spread out between the two sites; in the latter 
case, the EP coupling dominates and favors localizing the 
electron charge on only one of the two sites. The lattice 
spontaneously distorts so as to set up an attractive EP 
"potential well" which binds and localizes the electron. 
The electronic binding energy thus gained in turn sta- 
bilizes the local lattice distortion. This self-localization 
mechanism is the essence of polaron formation. 

Localizing the electron on either one of the two sites 
is energetically equivalent due to the reflection symme- 
try ((1,2) (2,1)) of the underlying Hamiltonian. At 
the level of the 0th order adiabatic approximation, this 
symmetry is broken in the 2-fold degenerate 0th order 
groundstates u- = ±u^\ The existence of two degen- 
erate local minima in Wq- can thus be understood as 
a direct consequence of the symmetry breaking which 
accompanies the self-localization transition. In the 1st 
order adiabatic approximation, the lattice kinetic energy 
restores this symmetry by inducing tunneling processes 
between the two potential minima, thus giving rise to 
a non-degenerate groundstate in which the two degener- 
ate 0th order states are admixed with equal probability 
weight. 

From the above discussion it is clear that such tunnel- 
ing processes are accompanied by a transfer of electron 
charge between the two sites. These lattice tunneling 
processes, occurring within the EP-induced multiple- well 
Born-Oppenheimer potential, constitute the basic low- 
energy mechanism whereby self-localized electrons can 
move through the lattice. At higher temperatures, ther- 
mally activated over-the-barrier hopping will dominate 
the polaron transport ]T^ , p3t which, again, can be de- 
scribed as a purely lattice dynamical phenomenon. Thus, 
within the framework of the 1st order adiabatic approxi- 
mation, polaron formation and polaron dynamics is fun- 
damentally reduced to a problem of non-linear lattice dy- 
namics. 



We now turn to the 1st order potential correction 
Wi„(u„) (|29| ) in the two-site problem, shown for several 
values of t/Ep in Fig. |l](b). Since Wi{u), according to 
( ^ ) , is controlled by the w-gradient of the electron wave- 
function I 5'(m)), we should expect it to exhibit peaks 
wherever | ^'(m)) varies most rapidly with u. In the two- 
site problem, this occurs at m_ = where | ^(u)) changes 
its character from being predominantly localized on site 1 
to being localized on site 2, as discussed above. For large 
I u_. I, I 4'(m_)) approaches a constant, either | 
or I vl/(^2)^^ hence VFi_(u_) for | m_ |^ oo. In the 

polaron regime Ep > Ep , the primary effect of Wi is 
to enhance the tunneling barrier separating the two po- 
tential minima. In addition, Wi-{u-) will also tend to 
shift the two polaronic potential minima further apart, 
thus causing the tunneling barrier to become wider than 
in the 0th order potential Wq. Both of these Wi-effects 
tend to suppress the tunneling rate through the barrier. 
Even though Vl^i(O) may be small compared to the 0th 
order barrier height A^o (^5|), its effect on the polaron 
tunneling rates can be quantitatively of some importance, 
since tunneling rates are typically exponentially sensitive 
to changes in the tunneling barrier. 

In the delocalized regime Ep < E^ the primary ef- 
fect of Wi is to soften the harmonic restoring force con- 
stant = d^_W^{0) by an amount 

Xi_= _T4^i„(m_ =0) 
out- 

Thus, Wi also lowers the critical Ep for the on-set of 
polaron formation. However, in the large-M limit where 
the adiabatic approximation is valid, that is, for <^ t 
(see below), these corrections are smaller than the 0th 
order results (^, ^) by factors of order (fi/t)^ x (iJp/t)^. 
Provided t and Ep are of comparable magnitude and Q ^ 
t (see below), Wi does not qualitatively alter the basic 
structure of the lattice potential W in either coupling- 
strength regime. However, Wi can become qualitatively 
important in suppressing certain non-adiabatic processes 
in lattice systems, as will be discussed in the next section. 

IV. VALIDITY AND LIMITATIONS OF THE 
ADIABATIC APPROXIMATION 

The basic criterion for the validity of the adiabatic ap- 
proximation is that the longest time scale of the elec- 
tronic motion should be short compared to the shortest 
time scale of the lattice motion or, equivalently, the low- 
est electronic frequency scale should be large compared to 
the highest lattice frequency scale. In the two-site prob- 
lem, the lowest electronic frequency scale is the excitation 
energy between the electronic groundstate | '^{u)) and 
the 1st excited state which is at least 2t (at u_ ~ 0) or 
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larger. The highest lattice frequency scale is the phonon 
energy Q and hence we expect the adiabatic approxima- 
tion to work, provided that 



n<s:2t 



(41) 



In the polaron regime Ep > Ep"^\ the lattice (not the 
electron !) motion acquires an additional, low frequency 
scale, given by the polaronic tunneling splitting 2tp be- 
tween the ground- and 1st excited states in the double- 
well lattice potential W{u). This tunnehng energy scale 
is typically smaller than or, at most, comparable to the 
bare phonon energy scale fl, given the conditions where 
a polaronic double-well forms in the first place. Hence, 
the basic criterion (^l|) applies in the polaronic regime 
just as well as in the delocalized regime, regardless of the 
electron-phonon coupling strength. Criterion (^) ap- 
plies even in the strong-coupling regime (^) where 2tp 
becomes orders of magnitude smaller than fi. 

Although the foregoing point has been established for 
some 30 years now |l^,g^,Q, we wish to strongly re- 
emphasize it here because a great deal of confusion about 
this has been created in the more recent literature on the 
two-site problem, as for example in Ref. ||2^. The basic 
error in some of the more recent work is to regard the 
polaron tunneling splitting 2tp, rather than 2t, as the 
lowest relevant electronic energy scale. Doing so, one 
then arrives at the much too restrictive validity criterion 



n < 2tp 



(42) 



If correct, this would imply that the polaron regime 
can not be treated in the adiabatic approx- 
imation, since typically ip ^ f2 even under the most 
favorable conditions. In the strong-coupling regime (^) 
where tp <^ Q the adiabatic approximation should break 
down completely according to ( Ji^ ) . 

The fundamental misconception here is that 2tp is of 
course not the lowest electronic energy scale, but rather 
represents an energy scale of the lattice motion, as dis- 
cussed above. The relevant lowest electronic energy split- 
ting, between the electronic ground- and 1st excited 
states at fixed lattice coordinate u is at least 2t in the 
two-site model, regardless of whether Ep is small or large. 

To illustrate this point, we have generated exact nu- 
merical solutions of the two-site problem using the full 
Hamiltonian H without any approximation, and com- 
pared to solutions of the 1st order effective adiabatic ef- 
fective Hamiltonian i?Ad = ^^m + W , corresponding to 
the effective action ^Ad from (|2^). For both the exact 
and the adiabatic problems, we have used a sufficiently 
fine discretization of the m_ coordinate and a sufficiently 
large cut-off at large U- to ensure a numerical accuracy 
of better than 1% in the calculated energy splittings over 
the entire parameter range studied. In Fig. ^ we show 
the logarithm of the polaron tunneling splitting 2tp, that 
is, the excitation energy from the ground- to the 1st ex- 
cited states of the full electron-phonon system, as a func- 
tion of Ep/n for t = I and four different EP couplings, 



Ep = 2.5, 3, 4 and 8 which are well inside the polaronic 
regime {Ep > i;^"'*^). 

In addition to the exact solution, we show two differ- 
ent adiabatic solutions in Fig. ^, one obtained with the 
full adiabatic lattice potential W = Wo + W\ , the other 
using only the 0th order potential, W = Wq. These are 
being referred to in the following as the "full" and as the 
"simple" adiabatic solutions, respectively. As expected 
from the Holstcin-Lang-Firsov strong-coupling expansion 
[|l|,||j2| and from semi-classical (WKB) arguments, the 
tunneling splitting at fixed Ep and t decreases exponen- 
tially with 1/0, as indicated by a roughly linear depen- 
dence of ln(2ip) on 1/ri in Fig. ^. 

Remarkably, the full adiabatic result agrees with the 
exact solution to better than 14% over a parameter region 
0.15^ < 51 < 0.5t wherein 2tp varies by more than 9 or- 
ders of magnitude, including the regime where 2tp is or- 
ders of magnitude smaller than fi. The simple {W = Wq) 
adiabatic solution reproduces the qualitative features of 
the l/ri- and £'p-dependence of 2tp quite well, but the 
quantitative agreement is noticeably worse than for the 
full adiabatic solution. The agreement between the full 
adiabatic and the exact results is all the more convinc- 
ing in light of the fact that the tunneling splitting is 
"exponentially sensitive" to small errors or changes in 
the wavefunction inside the tunneling barrier. Thus, 
our comparison of the tunneling splittings constitutes a 
much more stringent test of the underlying approxima- 
tions than a comparison of, say, low-lying-state expecta- 
tion values or wavefunction amplitudes. Other exact nu- 
merical results for the two-site problem, such as reported 
e.g. in Ref. |^^, are generally in equally good agreement 
with the corresponding adiabatic solution, provided, that 
is, one exercises enough care to use the proper adiabatic 
wavefunctions | (/)), Eq. (^), in carrying out the compar- 
ison. 

As expected from (|4l]), the agreement between adia- 
batic and exact results deteriorates at high phonon fre- 
quencies when J7 becomes comparable to t. As a practi- 
cal matter, even for fl = 2^, the agreement is still quite 
acceptable. For applications to lattice systems, it is of 
interest to explore in some detail how the adiabatic ap- 
proximation actually breaks down as one enters into the 
"anti-adiabatic" regime 



t < ?t 



(43) 



As a limiting case, we consider the ionic limit i ^ 0, al- 
ready discussed above. Here, the Holstcin-Hubbard prob- 
lem can be trivially solved exactly . Obviously there 
cannot be any electron tunneling between the two sites 
and the polaron tunneling splitting 2tp vanishes. 

By contrast, in the simple adiabatic approximation 
W = Wo, Wo- approaches the double-parabohc poten- 
tial (|3^) for t 0, which has a tunneling barrier of finite 
height and width. The simple adiabatic approximation 
would thus predict a non- vanishing finite tunneling split- 
ting 2tp > even for i = 0, a clearly unphysical result. 
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If instead one uses the full adiabatic approximation, 
with W = Wo + Wi, the correct qualitative physical 
behavior of 2tp is restored by the Wi-term, shown in 
Fig. 0(b): According to Eq.(^ the VFi-peak height (at 
fixed Ep and H.) diverges as t^, while at the same time 
its peak width vanishes, but only linearly in t in the limit 
t —^ 0. It is then easy to show that the transmission am- 
plitude through the W^i -barrier vanishes, that is, the bar- 
rier becomes impenetrable in the limit t ^ which forces 
2tp ^ for t ^ 0. Thus, as far as the tunneling splitting 
2tp is concerned, the full adiabatic approximation repro- 
duces qualitatively the correct physical behavior even in 
the extreme anti-adiabatic regime. 

The actual failure of the full adiabatic approximation 
in the t — > 0-limit is a more subtle problem. It consists of 
the unphysical constraint being imposed on the dynam- 
ics of the u_-coordinate by the impenetrability of the 
Wi-barrier. For t — > , th e Wi-barrier forces the lattice 
wavefunction (j>{u) in ( p4|) to vanish identically either to 
the right (u_ > 0) or to the left (u_ < 0) of the bar- 
rier. Thus, the amplitude for propagation from an initial 
M_ < to a final m_ > (or reverse) vanishes in the full 
adiabatic approximation at i = 0. In the exact solution 
of the t = problem, this constraint does not exist; the 
lattice is free to propagate with some finite amplitude 
from M_ < to M_ > 0. In the exact t = solution, 
the lattice dynamics is governed either by the left or the 
right parabolic well, W'^^^^ or W^^^\ corresponding re- 
spectively to the left-localized or to the right-localized 
electron states, | vf^^^)) or | discussed in Sec. HI. 



The problem with the adiabatic approximation is that 
the t = electron groundstate | ^'(u)) exhibits a level 
crossing and thus changes discontinuously at U- = 0, 
as discussed above. The adiabatic approximation, by 
construction, excludes transitions between, say, the elec- 
tronic ground- and 1st excited eigenstates. But this is 
just what happens at U- = in the t 0-limit: If 
the lattice coordinate crosses U- — from the left, say, 
under exact time evolution, the electron remains in its 
localized state | '^Z^^^^), which is the groundstate only for 
U- < 0, but becomes the 1st excited state when u- > 0. 
The adiabatic approximation on the other hand forces 
the electron to remain in the groundstate which changes 
discontinuously at — 0, from | ^^^'^^) to | v]/(^2)^^ 

In the two-site problem, the foregoing impenetra- 
bility constraint causes only a small error, of order 
exp(— E'p/n), in the low-lying lattice eigenstates and 
energies if the lattice oscillator zero-point amplitude is 
small compared to the double- well separation \/2up , that 
is, if r2 <C Ep. However, the impenetrability constraint 
may introduce a qualitative failure of the adiabatic ap- 
proximation if applied to large systems N ^ oo and tun- 
neling processes which transfer a polaron in a single step 
over large distances, as we will now discuss. 

Let us consider for simplicity the case of the Holstein 
model for just a single electron in a large lattice with 
sufficiently strong Ep to form a polaron. Suppose the 
polaron is localized at some site ^, say, and we want to 



study the tunneling barrier for transferring the polaron 
in a single tunneling step to a distant site C = C + i-S-, 
with Irl 3> a where a is the lattice constant. 



Let = (wi' 



(0 



denote that lattice configu- 



ration which minimizes Wo{u) and localizes the polaron 
around the "centroid site" ^ G {1 . . .iV}. That is, juf^l 
and the corresponding electron charge density {n£)^^"> are 
maximal at ^ = ^ and die out exponentially at large dis- 
tances \i — ^1 from the centroid. Likewise, let w^''^ de- 
note the lattice configuration which localizes the polaron 
around site C- By lattice translational invariance 



,(C) 



,(«) 



(44) 



if C = ^ + r. Notice that polaron formation breaks the 
translational symmetry of the lattice in the 0th order 
adiabatic approximation. As a consequence, Wq exhibits 
N degenerate local minima, corresponding to the N dif- 
ferent, but translationally equivalent u^^^ configurations 
on an Ai"-site lattice with periodic boundary conditions. 
This is the lattice analogue to the breaking of reflection 
symmetry in the two-site problem. 

Let u'^^^)(s) denote the linear path segment in the iV- 
dimensional u-space connecting u'^^ to 

i.e., 



,(C) 



(45) 



with s G [— ^,-f:|]. In the following discussion, we con- 



sider (|45|) as a representative of low-action tunneling tra- 
jectories connecting m^^^ to u^''\ The s-coordinate can 
thus be regarded as the lattice analogue to the u_ tun- 
neling coordinate (E3) in the two-site problem. Note in 



particular that Wo(u^'^^' (s)) has local minima at s = — ^ 
and s = +^ which must, by continuity, be separated by 
(at least) one intervening maximum, i.e., by a tunneling 
barrier. The simplest scenario, normally borne out in 
the numerical calculations discussed below, is that there 
is only one barrier maximum, by symmetry located at 
s = 0. Thus, along w('»^^(s), Wq has qualitatively the 
same structure as Wo-(u_) described above for the two- 
site problem. 

The first crucial point to note here is that the width 
of this tunneling barrier, that is, the Euclidean distance 
from w'-^-' to li'''' in their A^-dimensional M-space, 



(46) 



is finite and bounded by an upper limit doo which is in- 
dependent of the spatial distance |C~'?I — kl- Note that 
doo is independent of ^ or ^ due to lattice translational in- 
variance. Thus two polaron configurations m'-^-' and u'^-* 
are never further apart from each other than doo in u- 
space, regardless of how far apart their centroid sites ^ 
and ^ are in real space. 

The second important point is that the height of the 
0th order (Wo) tunneling barrier along u*^^^)(s) is also 
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bounded independently of lattice distances |C~CI- To see 
this, note that the EP potential Cue acting on the elec- 
tron is attractive, i.e., Cui < 0, for any u-configuration 
along the path u(^'>)(s) between s = and s = 1. Hence, 
the contribution to Wo{u) from Hc+Hc-phiu) is bounded 
from above by the electron groundstate energy of the 
undistorted lattice. Also, by an argument analogous to 
(^), the elastic energy contribution Hk{u) is bounded 
from above by ^H^iu'^^^). Both of these upper bounds 
are independent of |C — CI- 

The foregoing considerations suggest that a manifold of 
tunneling trajectories exists, sufficiently close to u'-''^\s), 
which will all connect u^^^ to u^^^ through a Wo-barrier 
whose height and width is bounded by upper limits inde- 
pendent of IC ~ CI- Within the simple adiabatic approx- 
imation, W = Wo, one thus arrives at the unphysical 
result that the polaron can tunnel in a single ("instan- 
ton") tunneling step from any site C to any site C in the 
lattice with a tunneling matrix element ip(C ~ C) which 
does not go to zero for |C — ^| ^ cx3, but rather 



lim |ip(C-C)l =ipoo > 



(47) 



The foregoing argument can be made formally more rig- 
orous, by employing instanton methods similar to those 
described in the next section for short-distance tunneling 
processes. We will not engage in that exercise here. Suf- 
fice it to say that the simple adiabatic result (|4^) for the 
lattice is analogous to the above described two-site result 
in the t — limit: the simple adiabatic approximation 
allows tunneling solely on the basis of the Wq electronic 
groundstate energy barrier, regardless of whether there is 
actually any electronic wavefunction overlap between the 
initial and final it-configurations of the tunneling process. 

To account for wavefunction overlap effects in long- 
distance tunneling processes, the Wi-term (22) has to 
be included in the total potential W = Wq + W\. Let 
us consider the evolution of the electronic groundstate 
wavefunction I '^[u)) along the linear tunneling trajec- 
tory (s) (BS) between two centroid sites C and C with 
IC — CI ^ ^p('")- Here Ipiu) denotes the exponential lo- 
calization length of I ^'(m)) for local lattice distortions 
comparable to u^^^ . As a simplest scenario, let us as- 
sume that the wavefunction | ^'(u)) remains localized for 
all u along u^''^)(s). This situation will be realized at EP 
coupling strengths -Ep which are sufficiently large com- 
pared to E^^^^^ . The electronic groundstate problem can 
then be qualitatively described as follows: 

The EP potential Cuf^^is), acting on the electron at 
sites £, consists of two localized wells, C(i — s)?/^^'' and 



C(| -I- s^u'P , the former centered around site I — the 
latter around t = C,. As s is varied from — ^ to -t-i, 
the EP well at C becomes shallower and the EP well at ^ 
deepens. At s = 0, the two wells become degenerate. As- 
suming large real-space tunneling distances |^ — C|, the 
electron wavefunction overlap between these two wells 
is exponentially small. Hence, the electron groundstate 



wavefunction | ^(u'^''^)(s))) will remain localized at site C 
for most s < until s gets very close to s = 0. Within a 
very small interval around s = 0, | {u^'^^\s))) will then 
switch over from being localized around ^ to being lo- 
calized around C- In that narrow s-region around s = 0, 
the electron wavefunction consists of the superposition of 
two almost non-overlapping localized parts, one centered 
around ^, the other around Q. Since | ^{u)) changes very 
rapidly as a function of u near w^^f (0), Wi{u) will exhibit 
a sharp peak along it^^^^(s) which increases the tunnel- 
ing barrier at s = and hence suppresses the tunneling 
amplitude. 

Formally, this problem can be treated by a tight- 
binding ansatz for the electron groundstate wavefunction: 
I ^'(it^''^-' (s))) near s = is approximated by a superpo- 
sition of I *(u(«V2)) and | ^'(u'^V2)), i-e-, by the single- 
well groundstates of the two EP wells ^Cu"P and ^Cu'p, 
discussed above. As s is varied near s — Q, the response 
of I \E'(m*^^^)(s))) to the changing EP well depths is then 
governed by the effective electronic hybridization overlap 



ieft(C-C)= (*(^^«V2) I i?c I *(w(^V2)) 

,IC-Ch 



f exp 



(-^) 

' 2 



(48) 



where Ip ^ = ip{u^^^ /2) is the localization length of 

I ^(u'-^V2))- Within the tight-binding ansatz, the prob- 
lem then becomes analogous to the two-site problem 
in the i — > limit, with the tight binding-basis states 
I "^{u^^^ /2)) and | vI/(-u(C)/2)) replacing the two-site ba- 
sis states I ^f^^^)) and | ^-t^^)^^ respectively. Wi{u'^'^^\s)) 
exhibits a sharply peaked barrier at s = 0, analogous to 
the t ^ Q limit of the two-site problem. The Wi -barrier 
will be roughly of the form given by Eq. (p9[), with u_ 
replaced by u-(s) = c'(C;C) ^ s and with t replaced by 
^cff(C ~ 0- Thus, along with tcff(C ~ C)j the transmis- 
sion amplitude through the Wi -barrier and the effective 
polaron tunneling matrix element tp(C ~ C) will decrease 
exponentially with the tunneling distance |C — CI? analo- 
gous to the t — > 0-limit in the two-site problem. 

The long-distance polaron tunneling processes are in 
the anti-adiabatic regime, since the relevant effective elec- 
tronic hybridization overlap matrix elements teff(C ~ C) 
become small compared to the phonon energy at large 
tunneling distances |C — CI on large lattice sizes N . The 
Wi potential ensures, at least qualitatively, that the ef- 
fective polaron tunneling matrix elements ip(C — C) are 
properly suppressed to zero at large tunneling distances. 
Hence, the full adiabatic approximation W — Wq + W\ 
restores the correct long-distance behavior, as far as the 
polaron tunneling matrix element is concerned. 

However, just as in the anti-adiabatic limit of the two- 
site problem, the Wi-term also imposes an unphysical 
constraint on the lattice coordinates. In the present case, 
involving long-distance tunneling on a lattice, this con- 
straint acts to couple the lattice displacement coordinates 
at arbitrarily large distances |C ~ CL thereby introducing 
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unphysical infinite-range interactions between the lattice 
coordinates. 

Thus, in long-distance tunneling processes, the pre- 
conditions for the adiabatic approximation break down. 
However, from the foregoing discussion it is clear that 
the effective action for the corresponding paths increases 
exponentially and that the corresponding tunneling ma- 
trix element dies out exponentially with the tunneling 
distance. The simplest way of dealing with such long- 
distance tunneling processes is therefore to altogether 
neglect the corresponding tunneling paths in the path 
integral. This is what we will do in the following anal- 
ysis. As far as the polaron tunneling dynamics is con- 
cerned, the short-distance processes will be dominant. 
The relevant effective electronic matrix elements icff for 
short-distance processes are of order of the 1st neighbor 
t which is normally larger than or at least comparable to 
the phonon energy scale in typical solid state situations. 
We can therefore use the adiabatic approximation to ac- 
curately estimate the effective action for short-distance 
tunneling paths. And it is only in this limited sense that 
the adiabatic approximation willhe used in the following. 



V. INSTANTONS AND EFFECTIVE 
HAMILTONIAN 

The problem of polaron formation in the 2D Holstein- 
tJ and Holstein-Hubbard models has already been stud- 
ied extensively. JTo|- p^ , p^ In the nearly ^ -filled band 
regime, the dopant induced hole carriers in the AF spin 
background can form polarons with much less EP cou- 
pling strength than is required for a single electron in an 
empty band. Thus Ep for forming a single hole po- 
laron in the ^-filled system is reduced by a factor of about 
4—5, compared to a single electron polaron formation 

fcrit) 

in the empty band system. This reduction in Ep has 
been explained in terms of the hole mass enhancement 
and self-localization effect in the AF spin background of 
the nearly i -filled Hubbard system [Q. The basic idea 
here is that the coupling to the AF spin background al- 
ready provides some form of self-localization of the hole 
relative to a self-induced local distortion of the AF spin 
correlations ||l|,|^. This spin polaron effect is mani- 
fested in the strongly reduced hole quasiparticle band- 
width, from 8t in the non-interacting system to ~ 2 J in 
Hubbard or tJ systems near half-filling. In the presence 
of EP coupling, this electronic bandwidth reduction per- 
mits the hole quasi-particle to become self-trapped by a 
much weaker EP potential well; hence the reduction in 

e: 



(crit) 
P 

p(crit) 



The fact that the polaron formation threshold 

Ep"""' > remains non-zero even in the strongly corre- 
lated systems is dictated by the so-called small-polaron 
dichotomy |2|], as discussed further in Sec. IX. 
For a multi-hole system containing 



P^N-J2n, (49) 



doping-induced holes on an iV-site lattice, there are (p) 
possible configurations for accommodating the P local- 
ized holes on the N available sites. The lattice poten- 
tial W{u) is therefore expected to have up to (p) nearly 
degenerate local minima, denoted by in the follow- 
ing, corresponding to the (p) different centroid configu- 
rations ^ = {^1, . . . ,^p) Here, = (6,x6,y) denotes 
lattice (centroid) site occupied by the i-th hole. As noted 
above, each of these local-minimum configurations breaks 
the translational symmetry of the lattice at the level of 
the 0th order adiabatic approximation. The symmetry 
is restored in the 1st order adiabatic approximation by 
polaron tunneling processes between the different u^. 

At EP coupling strengths Ep larger than, but suffi- 
ciently close to Ep''^^^ , it is possible that some of the (^) 
centroid configurations ^ do not have corresponding sta- 
ble local-minimum configurations in W{u). This may 
happen, for example, in a two- hole system {P ~ 2), if one 
tries to accommodate the two polarons at 1st neighbor 
sites, ^1 and ^2, in the presence of a 1st neighbor Coulomb 
repulsion Vc- At sufficiently strong Vc, the correspond- 
ing local minimum = u(^i'?2) becomes locally unstable, 
which is signaled by the smallest eigenvalue of the restor- 
ing force matrix d'^W/du^l^i becoming negative. In the 
following, we will not consider such situations, but rather 
restrict ourselves to parameter regions where all the local 
minimum configurations are stable. 

To establish the basic structure of the effective polaron 
tunneling dynamics, we treat the path integral for the 
effective action SAd ( p3| ) or its equivalent Hamiltonian 
^Ad by a lattice dynamical many-body tight-binding ap- 
proach. The basic idea behind this approach is that 
an effective polaron tunneling Hamiltonian Hp can be 
defined which operates in a "low-energy" sub-space of 
nearly orthogonal tight- binding basis states | 0^), la- 
beled by the localized polaron centroid configurations 
^. Each such basis state represents a lattice wavefunc- 
tion 4>^{u) which is assumed to be localized in u-space 
around the corresponding local potential minimum con- 
figuration u^. For example, <j>^{u) could be chosen as the 
vibrational ("phonon") groundstate obtained in a har- 
monic approximation by expanding W{u) to quadratic 
order around u^. By restricting the lattice Hilbert space 
to such a set of basis states all vibrational excited 
states around the polaronic local minima are neglected. 
Thus, formally, our approach can be regarded as a tight- 
binding approximation, formulated for the quantum dy- 
namics of the multiple- well lattice potential W{u) in the 
A^-dimensional lattice configuration (u-) space. 

In the simplest tight-binding approach one would then 
simply estimate the matrix elements of Hp by projecting 
the adiabatic lattice Hamiltonian H^d onto the corre- 
sponding tight-binding low-energy sub-space spanned by 
all ^j. In such a 1st order projection approach, one ne- 
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gleets all effects arising from virtual excitations out of 
the low-energy sub-space. 

Since tunneling matrix elements are exponentially sen- 
sitive to small corrections in, for example, the tunneling 
barriers, such a tight-binding projection could cause se- 
vere quantitative errors in the estimation of the magni- 
tude of tunneling matrix elements. Also, as a practi- 
cal matter, the accurate evaluation of Hamiltonian ma- 
trix elements with basis functions defined on the A^- 
dimensional it-space can become quite difficult. Lastly, in 
the 1st order projection approach, it is difficult to include 
the tunneling Berry phases into Hp. 

To avoid the foregoing difficulties, we have adopted 
an approach which is based on a direct mapping of 
imaginary-time tunneling paths, rather than a mapping 
of Hamiltonian matrix elements. Formally, this is accom- 
plished by the path integral instanton method |2^-|3^ . In 
this method, an "instantaneous" polaron hopping pro- 
cess ^ — > C induced by Hp between two polaron centroid 
configurations ^ = (^i . . . ^p) and C = (Ci . . . ^p) is identi- 
fied with the (restricted) path sum of instanton tunneling 
paths connecting to in it-space. The effective ac- 
tion Sp of the instantaneous hopping paths, so obtained, 
can then be immediately translated into matrix elements 
of the effective tunneling Hamiltonian Hp. Since only 
tunneling paths, but no basis states enter into the map- 
ping, the results do not depend on any particular choice 
of tight-binding basis states (j)^. 

As a specific starting point, we consider the trace of 
the resolvent operator at complex energies E 



which are almost always close to one of the centroid 
configurations, occasionally transfer from one to another 
centroid configuration by an almost instantaneous po- 
laron hopping process, and finally return to the initial u- 
configuration at imaginary time f3, in order to satisfy the 
closed path constraint. Important closed-path tunnel- 
ing processes for polaron states with P=l and 2 dopant- 
induced holes are shown in Fig. |[ Each black circle rep- 
resents an occupied polaron centroid site in the initial 
configuration ^ of the hopping process. Arrows indicate 
the hopping processes transferring the initial configura- 
tion ^ into the final configuration C. Thus, in u-space 
each arrow corresponds to a set of instanton-type tunnel- 
ing paths which connects the two respective minimum- 
W endpoint configurations and it^ and traverse the 
VF-barrier separating the two minima. Note that, 
as discussed above, via such tunneling paths, a hole po- 
laron can tunnel in a single step between 2nd-, 3rd-, etc. 
neighbor sites even if the original electron Hamiltonian 
(^ contains only a Ist-neighbor t. 

First, we consider the case of P=l. For the time be- 
ing, we take into account only the 2nd and 3rd neigh- 



bor processes denoted by amplitudes 4^"* ^^'^ ^i"' 
Fig. U(a). Single-polaron mier-sublattice processes are 
strongly suppressed by the AF spin correlations. ]3^ ] 
Hence, the Ist-neighbor amplitude i^^^ can be much 



,(3) 



smaller than or, at most, comparable to 4^"* and tf (to 
within 20-30%) in the case of P—1. Then, instanton path 
configurations are classified according to the numbers of 
intra-sublattice processes: Ux counts the number of tl 



A3) 



f°° (3) 

Tr (£' — H)^^ = dp e^^ Tr e~^^ , (50) processes to the right, the t\ processes to the left 

•^0 riy the tf^ processes to the upper, ruy the t^^"* processes 



written in the imaginary-time domain in path-integral 
form, 



Tr e 



Vu{t) 



~SAd["(T)] 



(51) 



u(/3)=«(0) 



(2) 

to the lower, n„ the t\ processes to the upper-right, m„ 
the t\ ' processes to the lower-left, the t\ ' processes 

(2) 

to the lower-right, and my the t\ processes to the upper- 
left neighbors. Path integration over the corresponding 
instanton paths gives 29 - 3l|] 



The trace operation in Eq. ( pOf ) leads to periodic bound- 
ary conditions on the imaginary time interval [0, /3] in 
Eq. (|5l|). These periodic boundary conditions in Eq. ( |5l| ) 
impose not only the closed path constraint u(t) — u(0), 
but also the condition that the initial and final elec- 
tron wavefunctions must be the same, including their 
phase factors. That is, for the electron wavefunctions 
I ^'(u(r))) entering into S'Ad via Eq. ([Tg|), the constraint 
(^'(u(/3)) I *(m(0))) = -f-1 must be imposed for aU paths 
u(t) integrated over in Eq. (|l|). The latter require- 
ment ensures that the Berry phase contribution to S'Ad in 
Eq. ( ^ ) is uniquely defined for every closed path w(t), in- 
dependent of the choice of phase for each individual elec- 
tronic wavefunction | \l/(u(r))) along such a path. Quan- 
tized eigenenergies can be found from ( |5C| ) by searching 
for the poles of the trace of the resolvent operator on the 
real i?-axis. 

The main contributions to the low-energy part of ( |50| ) 
arise from instanton path configurations, i.e., w-paths 



Tr e-f'^ = 



E 



dpx 



27r 



dp. 



y_ ^ipy (2ny — 2my —niu —n^+m^ ) 



27r 



dpx dpy 



271 2n 



-/3{m^(w('""'1)) + 24')[cos(p, +py) 
+ cos{px~Py)] +2t^i\cos{2px)+cos{2py)]j . (52) 



The effective hopping matrix elements are obtained 
as 
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t 



5R 



(")_ 



(53) 



the absolute minimum lattice potential en- 
ergy obtained at a minimum- configuration = 
for p = 1. Factorial factors such as rixl, etc. are 
introduced to account for identical species of instantons. 
The Px and Py-integrals are introduced to enforce the 
imaginary-time periodic boundary condition. The quan- 
tity SR'-p^ is the single-instanton contribution to the real 
part of the action for the path segment of the correspond 



function in the path integral takes the form of the trace 
over a 2 X 2 matrix exponential, namely 



Tr e 



-I3H 



dPx dpy 

ir exp 

27r 27r ^ 



-PW{u 



(min,2) \ 



1 
1 



~/3 



2*2^^ COSPx 

MSf^ cos ^ cos ^ 

(mill, 2) ^ 



4t 



(2) 



cos ^ COS ^ 



2t, 



(3) 

2 COSPj, 



(54) 



, ,. Xv) , Ay) . where VF(u'™™' M denotes the absolute minimum lattice 

mg tunnelmg process tp and Vp is the correspondmg . ^ r n o u • j (min 21 — ff f 1 



Berry phase contribution. The assignment of a unique 
Berry phase factor e^**J= to each such open path segment 
requires more detailed symmetry considerations and will 



potential energy for P = 2, obtained at u 



The quantity k'P 



^-th power of the fluctuation detcrmi- 



be postponed until Sec. |VI 
represents the 

nant for the instanton solution with the zero mode ex- 
cluded divided by that for the static solution at ^("""'i)^ 
and j'^p' is the Jacobian needed for a special treatment 
of the corresponding zero mode. They are defined as 
in Eqs. (10.13) and (10.14) of Ref. for the periodic 
potential problem. Substituting the result of the path 
integral ( p^ ) into the formula (^0|) , we obtain the disper- 
sion relation shown in the parenthesis of exp [— /3{. . .}] in 
Eq. (52). Note that the effective hopping matrix element 
is defined such that it is positive if the corresponding 
path segment carries a nontrivial (—1) Berry phase fac- 
tor: the sign convention of our polaron tunneling matrix 
elements is opposite to that used in the original elec- 
tron Hamiltonian (^. 

Next, we consider the case of P = 2. Since self- 
localization reduces substantially the polaron kinetic en- 
ergy scale, it is favorable for two polarons in an AF 
spin background to be bound in a pair: the binding en- 
ergy can easily be of the order of the effective polaron 
nearest- neighbor attraction, i.e., comparable to the AF 
spin exchange coupling J l33] in the Holstein-t J model. 
As a first approximation, we therefore restrict the path 
integration to include only 1st neighbor pair configura- 
tions u^^i'^^) and the instanton tunneling paths connect- 
ing them. Our numerical studies described below suggest 
that these 1st neighbor configurations represent the ah- 
solute minimum of W{u) for P = 2. Other, more distant 
pair configurations with |^i — ^2! > 1 are either repre- 
sented by local VF- minima it(?i?2) higher energy or they 
don't form local minima in W{u) at all. We are thus 
limiting ourselves, for now, to the tunneling processes 
and 

w-configurations as shown in Fig. §(b). 

The technique used above for P = 1 can be generalized 
in a straightforward manner to the present case P = 2. 
Here, in addition to the lattice translational degeneracy 
of the minimum-M^ u-configurations, the P = 2 system 
exhibits 2-fold internal degeneracy, corresponding to the 
two possible orientations of the 1st neighbor polaron pair, 
along either the x- or along the y-axis. Because of this 2- 
fold internal degree of freedom, the instanton exponential 



and between the degenerate, absolute-minimum 



(min,2) ^ y(Ji52) 

with ^1 and ^2 denoting 1st neighbor centroid sites. The 
tunneling matrix elements ip'* are expressed analogous 
to Eq. (m) in terms of the action contributions, fiuctu- 
ation determinants, and Jacobians of the respective in- 
stanton path segments. The 2 low-lying eigenenergies of 
the polaron pair at total momentum p = (px,Py) are ob- 
tained by diagonalizing the 2x2 matrix in exp [— /?{. . .}] 
of Eq. (|5|). 

The generalization of the foregoing path integral ap- 
proach to P > 2 hole polaron states is in principle 
straightforward, but becomes practically difficult to im- 
plement with increasing polaron number P. Analogous 
to (|5^), the approach leads to a momentum integral 
over the trace of a matrix exponential where the ma- 
trix dimension refiects the number of (nearly) degenerate, 
translationally inequivalent polaron centroid configura- 
tions (^1 . . . £,p) included in the tunneling analysis. In the 
following, we restrict ourselves to the cases P =0, 1, and 
2 which will allow us to extract the effective 1-polaron 
tunneling and 2-polaron interaction matrix elements. 

The low-lying tunneling eigenenergies identified by the 
foregoing instanton path integral method (and their cor- 
responding eigenstates) can be equivalcntly represented 
in terms of an effective polaron tunneling Hamiltonian 
Hp where each polaron is represented as a spin-i fermion. 
Hp is thus defined to operate in an effective spin-^ 
fermion Hilbert space with the effective fermions occu- 
pying sites on the 2D square lattice of the original 
EP Hamiltonian. A P-polaron centroid configuration 
(fi ■ • ■ Cp) is mapped onto the corresponding state of P 
site-localized fermions with minimum possible total spin, 
i.e., with Stot = 5 (0) for odd (even) P. The latter map- 
ping condition reflects the fact that the absolute elec- 
tron groundstates | vE'(w)), numerically calculated on fi- 
nite clusters, exhibit minimum total spin quantum num- 
ber. Notice however that by representing the polaron 
as an effective spin-^ fermion, we are actually including 
low-energy spin excitations into the effective Hamiltonian 
description. In order to derive the effective polaron spin- 
spin interactions, our adiabatic path integral treatment 
can be straightforwardly generalized to include restricted 
electron groundstates in Hilbert space sectors of higher 
total spin quantum numbers Stot > 1- In this manner, 
the spin-i fermion representation can be extended well 
beyond the scope of our original adiabatic approximation 
which retains only the (minimum-spin) absolute electron 
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groundstate | ^{u)). In the following analysis, we limit 
ourselves to the absolute groundstate only. Hence we 
are only studying the total spin-singlet pair state in the 
P = 2 case. Usin g ou r numerical Berry phase results, we 
will show in Sec. 



VI 



that each single polaron in such a 
singlet pair behaves indeed as a spin-i fermion. 

In generalizing the above 1st neighbor approach, it is 
also straightforward to include inter-sublattice hopping 
processes: the dimension of the matrix increases, the k- 
independent term is no longer proportional to the unit 



matrix, and t^p' (more precisely, tY', ^3 



,(1) Ala) 



and i''^) 



are defined as above. Then, the effective Hamiltonian 
describing the polaron tunneling dynamics and effective 
polaron-polaron interactions can be written in the form, 

Hp = ^ {Uj + ^ ^tijkridk) 

i^j,<y k 

X(l - ndj.-a)d]^d.i„{l - Udi-a) 

- ^ VpHdtndj ■ (55) 

f . 1 . 

Thus, dj^ creates a spin- 2 -fermion polaron with spin a 

at site j, Udj = J2a djadja = 0,1 and P = J^j ^dj = 1, 2. 
The hopping term is to include, appropriately, the am- 
plitudes (t,j + Efc Ati.fendfc) = t^i\ t^i\ 4'"\ 



tj^'', or t2^' [with appropriate sign according to the cor- 
responding Berry phase factor] for i j tunneling pro- 
cesses shown in Fig. |[ Note here that the sign convention 
for the polaron tunneling amplitudes tp "* in ( |55| ) is oppo- 
site to that used in the underlying electron Hamiltonians 

(ID- 

The 1st neighbor attraction Vp in ( [5 51 ) is estimated as 

Vp = 2W^(u("""'i') - M^(m("""'2)) - M^(m(™"^o)) , (56) 

where Ty(M(™"^-P)) is the (absolute) minimum potential 
energy W{u) for the u-configurations = u^.ii---^p) 

which minimize W{u) for P holes. For P~2, our numer- 
ical calculations suggest that the absolute minimum- T4^ 
u-configuration does indeed correspond to the 1st neigh- 
bor pair. For purposes of estimating Vp numerically on 
small model clusters, we have minimized VFq instead of 
W = Wq + Wi, thus neglecting the effect of Wi on 
y(min,p)_ the physical parameter regime of interest, 
VL <^ t,Ep, these W^i -corrections to indeed 
small, of order {n/t)^. The full potential W ^Wo + Wi 
was used to calculate VF(u'^™™'^^). 

To obtain order of magnitude estimates for tp\ 
we have used both the dilute-instanton-gas approach 
p9|-^, as explained above, and a constrained lattice 
dynamics approach ||11| w hich is more straightforward 
and adopted in Sec. [VII^ . The two approaches have 
given similar results. In the latter approach, the lat- 
tice Schrodinger equation corresponding to SAd is solved 
exactly for u constrained to the linear tunneling path 



,(3) 



u'^^^)(s) which is defined analogous to Eq. ( ^ ) and con- 
nects the two energetically degenerate, minimum- ly po- 
laron end-point it-configurations and of the respec- 
tive hop. The hopping matrix element \t^p^\ is then esti- 
mated as one half of the ground-to- Ist-excited state en- 
ergy splitting. 



VI. SYMMETRY OPERATIONS AND BERRY 
PHASES 

Before going into numerical estimations of effective 
model parameters, we need to settle the quasiparticle 
statistics and the signs of effective polaron hopping pro- 
cesses by calculating Berry phase factors. To calcu- 
late exp(— z0[u(t)]) for tunneling paths u{t) shown in 
Fig. 1^, we discretize r with at least 5 r-points per lin- 
ear path segment and obtain | ^'(u(r))) of the Holstein- 
tJ model by the Lanczos exact diagonalization method 
on an iV=4x4 lattice with periodic boundary conditions. 
The electron Hilbert space is restricted to the sector of 
minimum total spin {S=0, 1/2, for P=0, 1, 2, re- 
spectively), which comprises the absolute ground state 
I ^'(u)) for u-configurations near the local M^-minima. 
The results for all paths in Fig. are summarized by 



i{t)] = tt 



,(2) 



.(3) 



(57) 



where to'-'') is the number of i^-th neighbor hops with 
for P—2 denotes the number of 1st- 



=2, 3, and TOj^^ 



neighbor hops indicated by the dashed bonds shown in 
Fig. ||(a) by the first polaron in close proximity to the 
second, static polaron, indicated as a black circle. The 
effect of the TOj^'-term can be illustrated, for example, 
by comparing the Berry phase factors e~*^ of the trian- 
gular paths (a) (A) and (b)(B) shown in Fig. ^. In both 
paths, a single polaron is taken around the triangle in 3 
steps, consisting of two 1st neighbor and one 2nd neigh- 
bor transfer. For the one-polaron case, (a)(A), the phase 
factor is (—1), for the two-polaron case (b)(B) it is (-1-1). 
Thus, the close proximity of the second, static polaron 
in (b)(B) has altered the Berry phase of the first polaron 
tunneling around a closed path. 

The origin of the TO2^-'-term can be traced back to the 
internal symmetry of | ^>{u)): For the local minimum-W 
u-configurations of 2nd- and 3rd-neighbor polaron pairs, 
I ^'(u)) is odd under reflection along the dashed lines 
shown in Fig. ^(b), i.e., along the pair axis for the 2nd- 
neighbor pair and perpendicular to the pair axis for the 
3rd-neighbor pair. Suppose, for example, that the first 
polaron hops from (1,0) to (1,1) in a first step and from 
(1,1) to (0,1) in a second step with the second polaron 
staying fixed at (0,0). These are the first two steps of 
path (b)(B) in Fig. ^. Note that the two steps generate 
the same final centroid configuration as a refiection along 
the dashed (0,0)-(l,l) line, shown in Fig. |5|(b). Because 
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of this odd "internal" parity of | ^'(w)) for the intermedi- 
ate (2nd neighbor polaron pair) configuration, one of the 
two 1st neighbor hops must contribute an additional fac- 
tor (—1). Assigning this (—1) phase factor to one of the 
two 1st neighbor steps in path (b)(B) of Fig. || is to some 
extent arbitrary. The pattern of dashed-line and full-line 
bonds surrounding the static polaron in Fig. ^(a) rep- 
resents one possible assignment which is consistent with 
all the closed-path Berry phase results in Fig. ^(b). As 
a consequence of its odd internal parity, the 2nd neigh- 
bor polaron pair configuration is actually allowed to con- 
tribute with finite amplitude to polaron pair wave func- 
tions of (1^2 _y2 symmetry, in spite of the fact that the 
2nd neighbor pair axis points along the nodal axis of the 
dx2-y2 pair wavefunction. 

The m^^^- and m'^'^^-terms can be regarded as due to 
strong antiferromagnetic correlation. Suppose a polaron 
is initially located at (0,0) and hops to (2,0), (1,1), and 
then back to (0,0) along the path (a)(B) in Fig. |. The 
electron initially located at (2,0) hops to (0,0) and then to 
(1,1), while the electron initially located at (1,1) hops to 
(2,0). Thus, if one approximates the AF spin background 
by a Neel state, two electrons of like spin are exchanged. 
This produces a fermionic (— l)-factor. More generally, 
when a closed path consists of an odd number of 2nd- or 
3rd-neighbor hopping processes, an even number of elec- 
trons within a sublattice are cyclically permuted, pro- 
ducing the (—1) factor within the Neel approximation to 
the AF spin background. In order for this to occur, the 
AF spin correlation has to be strong, but it need not 
be long-ranged. For P=l, the Berry phase rule can be 
completely explained in this way. 

For both P=l and 2, 6'[m(t)] is given by a sum of 
independent single-polaron hopping contributions and 
exp(— i0[u(r)]) does not depend on whether or not the 
two polarons are being adiabatically exchanged in a given 
path [Fig. |(b)Lil Thus, for example, paths (b)(B) and 
(b)(C) in Fig. ^ contain the same 1st and 2nd neighbor 
hops and they have the same Berry phase. The two paths 
differ only in that (b)(C) exchanges the two polarons, 
whereas (b)(B) does not. Since the pair is a total spin 
singlet, this implies that each single polaron in the pair 
behaves effectively as a spin-1/2 fermion or as a spin-0 
boson. Only the spin-1/2 fermion representation is con- 
sistent with the half-odd-integer total spin in odd-P sys- 
tems and, as discussed in the previous sections, it is the 
one we have adopted. Equation ( p7| ) rules out the pos- 
sibility of representing dopant-induced hole polarons as 
spin-0 fermions or as spin-1/2 bosons. 

To settle the signs of effective polaron hopping pro- 
cesses, we need to define Berry phase factors for the 
corresponding single-hop open path segments. Let the 
initial u-configuration of such a single-hop path segment 
be denoted by m*-^-* and the final u-configuration by u^'^K 
The assignment of a Berry phase to such a path seg- 
ment can be made unique by fixing the phase of the 
corresponding wavefunction | \E'(u'^^))) relative to that 



of I ^'(m^^'')) in some unique manner. Given w^^^ and 
I *(w(^^)), let I ^'(■■<=f)(M(^')) denote such a final-state 
reference wavefunction. Also, let | vf^'^'^) («('»))) denote 
that groundstate wavefunction | ^1/(^(0)) which one ob- 
tains by adiabatically evolving | ^(u)) along the tun- 
neling path segment, without discontinuity in phase, be- 
ginning with I ^'(m^^^)). The Berry phase of the path 
segment is then defined as the phase difference between 
I *(^d)(y(c))^ and I ^^^f^w^^')): that is, as the phase of 
the wavefunction overlap {<b^'''^\u^'^^) \ (m^^')). If, 

for example, m^^^ and m^'"-' are related by a lattice symme- 
try operation, we can choose | 'J'^^^) ('"'■'''')) as the ground- 
state wavefunction | ^'(u^^))) generated by applying to 
I ^(u'^^))) the symmetry operation which transforms m^^-* 
into u^''\ If M*^^-* and u'^'') are related by several differ- 
ent symmetry operations giving different | ^'-'''^\u'-<^)) , 
we need to specify the reference | ^(u''"-')), i.e., which 
symmetry operation is chosen to generate the reference 
I ^(u'^''))) from I 4'(m^^^)). There does not always ex- 
ist such a symmetry operation to relate | ^'(■u^^))) and 
I 4'(m^^^)), e.g., the 2nd- or 3rd-neighbor pair and the 
Ist-neighbor one. Then, we can arbitrarily choose the 
phase of | Vt^™^^ (m^^')). The Berry phase factor for the 
corresponding path is also arbitrary. Figure ||(a) is an 
example. If we chose the different phase (i.e., the nega- 
tive) of I *("f)(u('^))) for the 2nd (3rd)-neighbor pair, the 

signs of all the 4^"^ (^2^^^) processes would be reversed. 

For ^=1, we first fix, arbitrarily, the phase of | 
*(u'*))) for the centroid configuration ^=((0,0)). All 
the I ^'(''''^'(w'''^)) are then uniquely defined by either 
translation or rotation operations. The Berry phase fac- 
tors for P=l are summarized in Fig. ^(a). Here, the 
translation (a;, y) (x + a,y + b) is denoted by T(a, b), 
and the rotation (x,j/) — s- (xcoscj) — ?/ sin 0, a; sin -|- 
ycos(j)) is denoted by i? (</>). The left-hand side shows 
I vl/(ad)^y(C)^^^ xhe right-hand side shows the possi- 
ble choices of | \E'^''®^-'(u^^^)), generated from the same 
I 4'(m*^^^)) by the appropriate lattice symmetry opera- 
tions. Note that, for 2nd and 3rd neighbor hops, both 
translation and rotation operations generate the same 
I (m^^))). The first three lines of Fig. |(a) give 

eii-p[-ie[^^]^+l, exp[-i6'f^]=-l, exp[-i6'f^]=-l, and 
thus 

4'^ < , 4'^ > , 4" > ■ (58) 

Obviously, these results are consistent with the relation 

For P=2, we first fix the phase of | 4'(m^^^)) for the 
centroid configuration ^ = (^1, ^2) = ((0,0),(1,0)). Rotat- 
ing this I *(ii)) by angle 7r/2, we define the | vl/(''°f) (^(0)) 
for C = (Ci,C2) = ((0,0),(0,l)). The I vl/(-f)(w«))) for 
the other 1st neighbor pair configurations ^ are defined 
by applying translation operations to either | 5'(it'^^^)) 
or to its rotated version R{tt/2) \ '^(yS^'^)). The re- 
sulting Berry phase factors for P—2 are summarized in 
Fig. @(b). The first two lines imply that exp[— i6'2^-']=— 1, 
exp[— ^i^d thus 
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i?^ > , t;,''^ > 



(3) 



(59) 



The last two lines in Fig. ^(b) imply, by comparison to 
the first two lines, that available rotations would generate 
the same | (u^^^)) as the translations. The 1st- 

neighbor hops fj^"' 



and 1^2^"^ are positive or negative for 



the processes indicated by the dashed and full bonds, 
respectively, of Fig. ^(a), as discussed above. 



VII. BERRY PHASES AND QUANTUM 
NUMBERS 



Using the effective Hamiltonian (|55|) with parameters 

t^^\ 4^°\ 4^^\ 4^^ (with signs determined 

above), and ^p, we can now calculate the low-energy 
eigenstates for the P = 1 and P = 2 polaron systems. 

In the case P=l, the dispersion relation from Hp is 
given by 



ei(p) = 2t']^\cospx + cos py] + 2tY'[cos{px + py) 

+ cos{pa; - Py)] + 24^'' [cos{2px) + cos(2py)] . (60) 

As mentioned in Sec. on finite lattices, 4^"* is smaller 
than the 2nd and 3rd neighbor t's for P = 1. As the clus- 
ter size increases, the overlap between the two minimum- 

14^ wavcfunctions connected by the 4^' process, | ^'(u^^))) 
and I 'I'(m*^^^)) for a 1st neighbor bond (C,C)j becomes 
small. Then, the potential energy W{u) would develop a 
higher barrier for the 1st neighbor hop, due to Wi(u), so 

that 4^^ would become further smaller. Allowing for ar- 
bitrary values of tp' and tj^'' but 4^'' = 0, the one-polaron 
band-minimum is located at momentum p — (7r/2,7r/2) 



(2)r 



for I 4^^ |< 2tf\ at (7r,0) for | 4^^ |> 24^^ and 4^^ > 0, 
and at (0,0) and (tt, tt) for | 4^^ |> 24^^ and 4^^ < 0, 
as shown in Fig. |^(a). For the physically relevant signs 

(2) (3) 

implied by the Berry phase factors, t\ ,t\ > 0, the 
momentum of the one-polaron band-minimum is thus at 
(7r/2,7r/2) or (tt, 0) which lies on the Fermi surface of 
the noninteracting nearest-neighbor tight-binding band 
model at half filling. The one-polaron bandwidth is given 

by 



44^' + 8ti' 
84^) 



(3) 



for < <f ^ < 24^^ 
for < 24^^ < t^i^ 



(61) 



For the cluster geometries studied here {N = \/8 x \/8, 
N — vTO X -x/lO and iV = 4 x 4) and with only nearest- 
neighbor terms {t, J) included in the original Hamil- 
tonian (^, certain "accidental" symmetries exist which 

cause tp^ = tpK As a consequence, the band-minimum is 
at (7r/2,7r/2), and the eigenvalues of the inverse effective 
mass tensor at this point are 



(mri). = 84^' + 44^Ump) 



(62) 



(771^^)0 = 84' 



(3) 



At 



(2) 



At 



(2,3) 



(63) 



where the subscript r is for the (1,1) direction and is 
for the (1,-1) direction. If we include finite and negative 

4^'' (representing, e.g., the N —> 00 limit at finite, fixed 
hole density, rather than at fixed hole number P = 1), 
then the band-minimum is shifted by 4^' from {tt/2, n/2) 
to some point {p,p) with p < 7r/2 which would fall on 
the Fermi surface of the noninteracting band model at 
corresponding filling. 

For P—2, we first consider the tightly-bound pair limit, 
Vp S> I t'^'^ I , where we can approximate the polaron pair 
ground state by including only nearest-neighbor pair con- 



figurations, thus retaining only the and t]^' matrix 
elements of Hp. The pair dispersion relations are then 
given by 



,(3) 



e*(p)= — Vp + 4^^(cosp^ -I- cospy 



± 



t^^^ {cos Px — cos Py )'' + (4t2^^ cos ^ cos ^) 



(2) 



Py \2 



1/2 



(64) 



(2) (3) 

Allowing for arbitrary values of ^2 and 1 the pair 
wave function in the nearest-neighbor pair approximation 



> t^^'' has d^2_y2 -wave symmetry if t2 > 



(2) 



for I 42) 

and s-wave symmetry if ^2^'' < 0, and, in either case, 
total momentum p — (0,0), as shown in Fig. ^b), at 



the band-minimum. For | t 



(2) 



< t 



(3) 



the pair ground 



states are multiply degenerate: the horizontal pair (with 
pair axis parallel to the x axis) with total momentum 
p = {TT,Py) for arbitrary | py \< tt, and the vertical 
pair (with pair axis parallel to the y axis) with total 
momentum p = {px,Tr) for arbitrary \ p^ |< tt all have 
the same energy. The two-polaron bandwidth is given by 



Bo 



max €2 (p) — mine2 (p) — 4 
p p 



I -t 



(3) 



(65) 



If we take account of 1st-, 2nd-, and 3rd-neighbor pair 
configurations, including also the Ist-neighbor hopping 
terms, tj^""* and t'^''\ in second-order perturbation the- 
ory, the initially degenerate energy along p = (n^py), 



'iTT,Py) 



-Vp - 2t 



(3) 



is lowered by 



(^£2 (T^iPy) = - 



4(2t 



(i")2cos2ai+4")^) 



Vp + 2t 



(3) 



(66) 



(67) 



It is reasonable for the tj^^^-term to favor py = because 
the process of (a, ^2)=((0, 0), (1,0)) ^ ((0, 0), (1, 1)) and 
the process of ((0, 0), (1, 1)) ^ ((0, 1), (1, 1)), for exam- 
ple, are in phase, as shown in Fig. ^(a). The ground 
states are still doubly degenerate: the horizontal pair 
with p = (7r,0) and the vertical pair with p = (0,7r), 
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which would correspond to p^- and p^-wave though they 

are total-spin-sing lets, m For 4^\4^^ > implied by 
the Berry phase factors, we thus get either dj.2_y2- or 
P3.(j,)-pairing symmetry with total momentum p = (0, 0) 
or p = (tt, 0)[(0, tt)], respectively. 

The accidental symmetries for our finite cluster geome- 
tries (in the absence of longer-range terms in the original 
Hamiltonian (P studied here) lead to tg^'' = 4"^^ which 
is exactly on the d-p phase boundary where B2 vanishes 
due to a frustration effect. [|2| So, the energy e^(p) is 
independent of p. If we take account of 1st-, 2nd-, and 
3rd-neighbor pair configurations again, in second-order 
perturbation theory, the initially degenerate energy on 
the d-p phase boundary. 



■(P) 



,(,3)=-yp-24'''^ 



is lowered by 



(P) 



/(P) 



Vp + 21 



(2:3) 



(68) 



(69) 



with 



/(P 



1= 4t2^°^^(2 + COSP3. + cospy) 

_^^^(16)2 1 " COS Pa: COS Py 
^ 2 + COSpx + COSpy 



Then the ground state has d^^-y^ symmetry with p = 

(0,0) for \/2 I 4^"^ l>l 4^^^ I ^i^d P2;(y)-wave with 
p = (tt, 0) [(0,7r)] otherwise. For the iV=4x4 Holstein- 
tJ cluster with periodic boundary conditions, an acci- 



dental symmetry leads to | ^2 



(la) 



= 1 t 



(lb) 



and thus 

(la) 



d2;2_y2 -pairing symmetry. It is reasonable for the 
term to favor the dx2_y2-wecve state because the pro- 
cess of (^1,6) = ((0,0), (1,0)) -> ((0,0), (1,1)) and the 
process of ((0, 0), (1, 1)) ((0, 0), (0, 1)), for example, 
have opposite signs. Also, it is reasonable for the 4^''''" 
term to favor the pa;-wave state because the process of 
(^1,6) = ((0,0), (1,0)) ^ ((0,0), (2,0)) and the pro- 
cess of ((0, 0), (2, 0)) ^ ((1, 0), (2, 0)), for example, have 
opposite signs, as shown in Fig. ||(a). Once again we 
note that the 2nd-neighbor polaron pair configuration 
contributes to the polaron pair wave function of d,j,2_y2 
symmetry. 



VIII. EFFECTIVE HOPPING AND 
ATTRACTION 

We have seen how total momenta and internal sym- 
metries of few-hole-polaron states are determined by the 
signs and relative magnitudes of the effective polaron 
tunneling matrix elements. In this section, we show nu- 
merical estimates of them with effective polaron nearest- 
neighbor attraction and effective pair binding energy to 
see the energy scales of polaron dynamics. The relative 



energy scale of kinetic energy to interaction strength is 
controlled by the phonon frequency in the original Hamil- 
tonian (|^). It is noted again that we use a constrained 
lattice dynamics approach and exactly solve the lattice 
Schrodinger equation corresponding to the effective ac- 
tion (|2^ ) for the lattice-displacement configurations con- 
strained to the linear tunneling path of the respective 
hop, as described at the end of Sec. |^. The effective 
lattice potentials are based on Lanczos calculations on fi- 
nite clusters with periodic boundary conditions. The nu- 
merical results should be regarded as very rough order- 
of-magnitude estimates only. The nearest-neighbor at- 
traction Vp is calculated according to the formula ([5^). 
The pair binding energy A is estimated in the nearest- 
neighbor pair approximation, according to 



2ei(pi 



(min)x 



— / (miii)N 
£2 (P2 ) 



(71) 



where ei(p) and e2"(p) are defined in Eq. ( |60|) (with 

4^''=0) and Eq. (|64|), respectively, measured relative to 

the P = groundstate energy, and Pp"'"^ are the re- 
spective (thus different) momenta at the band minima 
discussed above. Note that the sign of A is so defined 
that A > signifies a net attraction, A < repulsion. 
Figure ||(a) shows the logarithm of the dominant 2nd- 



f 2") 

(70) and 3rd-neighbor hopping amplitudes tp and t^p for 
P=l, 2 in the Holstein-i J model on an = VSx a/S clus- 
ter. As expected in a polaronic system, ||l^ all tp-* are 
suppressed, roughly exponentially, with increasing Ep /ft 
and strongly reduced compared to the bare electronic 



,(3) 



p(crit) 



the tp^ can become 



t. However, for Ep near Ep 
comparable to the phonon energy scale fi. For P=2 
the proximity of the second, static polaron strongly en- 
hances the amplitudes tj^"* and 4^^ relative to 4^'' f^nd 

(3) 

t\ . It is worth noting that this effect occurs only in the 
presence of strong electron correlations where bipolaron 
formation is prevented by the strong on-site Coulomb re- 
pulsion. By contrast, this effect never occurs in ordinary 
polaronic systems with the electron-phonon interaction 
Ep larger than the local Coulomb repulsion. In the lat- 
ter case small bipolarons will form | |3^ which are much 
heavier than polarons. To generate the above-described 
delocalization (and hence mobility !) enhancement effect, 
it is essential to keep the two polarons spatially separated 
by strong enough on-site Coulomb effects. 

(2) (3) 

The accidental symmetries leading to tp = tp will 
be lifted on larger lattices and, more importantly, by in- 
clusion of longer-range couplings, such as 2nd- neighbor 
hopping t', Eq. (^0|), and extended Coulomb repulsion 
Vc, Eq. (prj), in the original EP Hamiltonian (jl]), as will 
be shown below. Due to the exponential dependence of 
the delocalization matrix elements ip on the lattice po- 
tential parameters, such additional couplings can sub- 
stantially affect the magnitudes of the tp'' parameters, 
without necessarily altering the Berry phase factors or 
the predominance of the 2nd- and 3rd-neighbor hopping 
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terms {tp' > tp ) and their two-polaron enhancement 

(4"^ > ii''^ ) ■ The Berry phase factors should be a robust 
feature of our model, since they reflect the topological 
properties of the relevant tunneling paths relative to cer- 
tain singular manifolds of the lattice action in u-space. 

Figure ||(b) shows the nearest-neighbor attraction Vp 
and the pair binding energy A, where the latter quantity 
is given by 



(2,3) 



8q 



(2:3) 



(72) 



for = t^p', using (|7l|). Since two self-localized nearest- 
neighbor holes mutually inhibit their delocalization, the 
t-term in the original Hamiltonian (|l|) gives a repulsive 
contribution to Vp : in the parameter range shown in the 
figure, Vp < 0.342J (0.316J) is substantially reduced 
compared to Vp(t = 0) = 1.00 J (0.926 J) oniV = V8x^/8 
(-\/T0 X -\/T0) sites in the i — > limit. Compared to the 
tJ model, Vp can be larger or smaller: self-localization 
reduces the effective polaron hopping processes, giving 
an attractive contribution, and it is more effective in the 
one-hole state than in the two-hole state, giving a re- 
pulsive contribution. The binding energy A is enhanced 

(2 3) 

by the two-polaron hopping amplitudes t2 , but it is 
smaller, in most of the parameter range shown in the 
figure, than Vp due to the restricted hopping processes 

(2 3) 

for the polaron pair and due to the non-negligible t\ ' - 
term for large fi. In a more realistic theory, the possible 
competition between polaron pairing and phase separa- 
tion would need to be considered for finite density of 
holes. 

In order to see a finite size effect, we have calculated 
the effective model parameters on iV = VTO x \/T0 sites 
(not shown) to compare with those on TV = \/8 X ^/8 
sites above. We find no qualitative difference between 
them. In the parameter range shown in the figure, the 
values of tp'' are different by a factor of two at most, 
but these values are rough order-of-magnitude estimates 
in any case. The values of Vp for N = \/T0 x \/T0 are 
smaller by a factor of 0.8-0.9. 

For the Holstein-Hubbard model, we find results 
(Fig. 1^) quite similar to those shown above. However, 
the values of Vp are only 30% of those in the Holstein-t J 
model, which are reminiscent of the fact that the hole 
binding energy is larger for the tJ model than for the 

(2 3) 

Hubbard model. Furthermore, the values of ^2 
smaller than those of the Holstein-t J model for 17 < 0.2t, 
and the values of 4^'^'' are larger by a factor of 1.6-3.7 in 
the parameter range shown in the figure. All these results 
make the pair binding energy smaller in the Holstein- 
Hubbard model. For large SI, the polaron pair becomes 
unbound, though our results are based on the adiabatic 
approximation and the nearest-neighbor pair approxima- 
tion so that they are less reliable for large O. 

We now turn to the effects of 2nd neighbor electron hy- 
bridization and long-range Coulomb couplings which lift 

(2) (31 

the accidental finite-cluster degeneracy, tp — tp' , and 



thus shift the system off the d-p phase boundary for P=2, 
already in the absence of 4^"''^''^-processes. The 2nd- 
neighbor electron hopping term in the original Hamilto- 
nian (|^) enhances the 2nd- neighbor hopping t^2\ lowers 

(3) 

the 3rd-neighbor one ^2 i and thus favors the d,j,2_y2- 
wave syrnmetry if t' is positive by the definition in Sec. ^ 
(Fig. 10), and the effects are opposite if t' is negative 
(Fig. 1^. Note that, in the noninteracting tight-binding 
model, the positive t' raises the energy of p—{'Kfi) state 
[thus the energy of p=(7r/2,7r/2) state is relatively low- 
ered] and makes the Fermi surface convex. A t'-term 
which helps the 2nd-neighbor electron hopping also helps 
the 2nd-neighbor polaron hopping. 

(3) 

The long-range repulsion term enhances t^ more than 
t^\ so that it favors the P2;(y)-wave symmetry (Fig. ^2]). 
This can be understood if we recall the second-order per- 
turbation theory with respect to t2^°'^^^/Vp. The Vq- 
term raises the energy of the intermediate 2nd-neighbor 
pair favoring the (ia,2„y2-wave symmetry, compared to 
that of the intermediate 3rd-neighbor pair favoring the 
P3,(j,)-wave symmetry. Note that Vc enhances both of the 

(2) (3) 

and t2 processes. This happens because the lattice 
distortion and thus the localization potential is weakened 
by Vq. If Vc is too strong, however, it may overcome the 
nearest-neighbor attraction and the polaron pairing will 
then be suppressed altogether. This will be discussed 
further in the next section. 



IX. POLARON LIQUIDS AND THE CUPRATE 
SUPERCONDUCTORS 

To the extent that the qualitative features of the above- 
discussed effective Hamiltonian (^) and the resulting 
tunneling and pairing dynamics remain intact at finite 
hole doping concentrations, the foregoing results have 
some potentially interesting consequences for the phys- 
ical properties of the polaron liquid formed at finite po- 
laron densities. In the present section, we will speculate 
on some of these properties and compare to experimental 
observations in the cuprate high-Tc superconductors. [T9| ] 

If the above-discussed polaron pair state remains stable 
and delocalized at finite hole doping, then formation of 
a superconducting polaron pair condensate can occur at 
low enough temperatures. The foregoing discussion has 
focussed primarily on the tightly-bound-pair limit where 
such a condensate would be formed via Bose condensa- 
tion of the pre-existing polaron pairs. However, the qual- 
itative Sl-dependences of the delocalization energies ip 
and of the pairing potential Vp, shown in Fig. ^ suggest 
that with increasing 17 (and fixed electronic parameters 
t, J and i?p), such a condensate may exhibit a crossover 
from tightly-bound-pair to a BCS-like, extended-pair be- 
havior: For small 12, the delocalization matrix elements 
t^2^ and resulting polaron pair bandwidth become small 
compared to the pairing potential Vp. Hence tightly 
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bound pairs will form, as described above, with a pair 
wavefunction extending only over 1-2 lattice constants. 
For large fl on the other hand, the polaron bandwidths 
(Bi and B2) can become comparable or larger than the 
pairing potential, thus leading to a BCS-like extended 
pair state, with a pair wavefunction extending over sev- 
eral/many lattice constants. 

In the tightly-bound-pair regime, the Bose condensa- 
tion Tc is controlled by the pair density Xpair and the pair 
bandwidth B2, that is, roughly 



pair J 



where B2 is the pair bandwidth (65) and 



1, , 1 

Cpair = T;il - W) = 



(73) 



(74) 



is the pair concentration, i.e., the half of the hole concen- 
tration X. i?2(^^) and hence Tc decreases with decreasing 

n. 

In the BCS-like extended-pair regime, Tc is controlled 
by the pair binding energy A which decreases with in- 
creasing delocalization energy and hence with increas- 
ing Q. As a consequence, there must exist, somewhere 
in the cross-over regime between the tightly-bound-pair 
and the BCS (extended-pair) limits, an optimal phonon 
frequency fio where the transition temperature Tc(il) is 
maximized, flo is roughly determined by the condition 



B2{no) ~ Vp 



(75) 



and the maximum possible Tc (as a function of f2 !), es- 
timated by extrapolation of ( [73|) from the tightly-bound- 
pair side, is of order 



T 



B2{^o) 



(76) 



where i?2(f^) is the polaron pair bandwidth, Eq. (|65|), as 
a function of phonon frequency fl. 

One crucial, experimentally observable difference be- 
tween the tightly-bound- and the extended-pair conden- 
sate is the relation between pair formation and supercon- 
ducting transition: In the tightly-bound-pair regime the 
pairs, and hence the pairing gap A in the excitation spec- 
trum, can be preformed. That is, the polaron pairs and 
the energy gap for pair breaking exist already at temper- 
atures T ^ A which could be well above Tc, provided 
that A ^ Tc- By contrast, in the extended-pair BCS- 
like regime we expect the pair formation to coincide with 
the superconducting transition, that is, the pairing gap 
should be observable only at temperatures T below Tc 
and should vanish at Tc- 

The existence of such an optimum phonon frequency 
implies that Tc exhibits a vanishing isotope exponent 
a when Q = fig- To show this, we note that the iso- 
topic mass dependence enters into the theory only via 
the phonon frequency fi, if the electron- phonon Hamil- 
tonian is parametrized, as in (^ and (^, in terms of 



Ep and fl, since electron-phonon potential constants (C) 
and harmonic restoring force constants (K), and hence 
Ep, are of purely electronic origin, i.e., do not depend on 
atomic/isotope masses. Using n oc Af-i/^ fr om p), we 
conclude that 



dlogTc 



dlogM 



^ I dlogTc 
2i 2(91ogf] 
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(77) 



which vanishes at the T^-maximum Q — flo. The notation 
... lei here means that the derivatives are to be taken with 
all purely electronic model parameters {t, U, J, Ep, etc.) 
held constant. The Tc-maximum at flo also implies that 
a is positive in the tightly-bound-pair regime flo > 
but negative in the extended-pair regime flo < The 
vanishing of a at = Qo does however not imply that 
a is generically a small number. Quite to the contrary, 
because of the strong fi-dependence of the polaron band- 
width parameters, we should expect a to attain quite 
substantial magnitudes, with |q;| ~ 0{1), as the system 
is tuned away from the optimal phonon frequency, i.e., 
when Q Qq. 

It is tempting to compare the foregoing features of a 
finite-density polaron liquid to the observed properties 
of the cupratcs. The doping-depcndence of the super- 
conducting and normal-state properties of the cuprates 
is, in some respects, very much reminiscent of a cross-over 
from tightly-bound-pair to BCS/extended-pair behavior: 
In the undcrdopcd cuprates, there is now a substantial 
body of evidence suggesting that the superconducting 
gap is pre-existing, in the form of a "pseudo-gap" , at 
temperatures well above Tc- With increasing hole 
doping concentration x, Tc approaches a maximum, while 
the pseudo-gap above Tc is gradually suppressed, and, in 
close proximity to the optimal doping concentration Xo, 
the pseudo-gap above Tc vanishes. Well inside the over- 
doped regime x > Xo, there is no detectable pseudo-gap 
and Tc rapidly decreases with increasing x. 

The isotope exponents a in the underdoped cuprates 
are typically quite large in magnitude, of order of the clas- 
sical BCS'-value aecs = 5 or larger. However, in con- 
trast to conventional BCS-type phonon-mediated super- 
conductors, a can be very sensitive to changes in doping 
and other system properties such as impurity concentra- 
tion and crystal structure. With increasing hole-doping 
concentration, the observed, usually positive oxygen iso- 
tope exponent a decreases and becomes very small, ty' 
ically < 0.05, at the optimal doping concentration Xo 
It is presently not clear whether a changes its sign for 
X > Xo- Negative a- values have been observed in cop- 
per isotope substitutions on less than optimally doped 
cuprate materials. 

In comparing these experimental results to the fore- 
going theoretical picture of a polaron liquid, it is im- 
portant to note that, experimentally, the Tc-maximum 
and the surmised cross-over from tightly-bound-pair to 
extended-pair BCS-like behavior is observed as a func- 
tion of doping concentration x, whereas, in our above 
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theoretical considerations, we have discussed the cross- 
over as a function of phonon frequency fl. To see how 
such a cross-over could arise in our model as a function 
of doping, we need to consider the doping dependence of 
the polaron delocalization matrix elements tp^ . 

As indicated in Figs. ^ |l2|, the polaron delocalization 
matrix elements, and hence the polaron pair bandwidth 
B2 are rapidly increasing functions of il. At finite dop- 
ing, these delocalization matrix elements will also become 
dependent on the hole doping concentration x = 1 ~ (n) 
by the following mechanism: As the polaron density in- 
creases, the localized wavefunctions of nearby holes will 
begin to overlap and the holes will begin to mutually 
screen out each other's tunneling barriers. This effect can 
be clearly seen in comparing the 1- and 2-hole results in 
Fig. ||. For P = 2, the mere proximity of the second, 
static polaron strongly enhances the tunneling matrix 
element of the first, moving polaron, hence t^'' > t'"-^^ 
for V = 2,3. Treated at the mean-field level, at finite 
polaron density, this tunneling enhancement effect will 
cause the (mean-field average) tunneling matrix elements 
to increase with the hole doping concentration. Thus the 
effective polaron pair bandwidth B2 = B2{^,x) becomes 
a strongly increasing function of of the hole doping con- 
centration X. 

According to the cross-over criterion ( [75| ) it may then 
be possible to drive the polaron liquid from the tightly- 
bound-pair into the extended-pair regime by changing ei- 
ther the phonon frequency f2 or the doping concentration 
X, if Vp is only weakly dependent on and x. Another 
way of stating the same result is to say that the optimal 
phonon frequency Qq ~ ^oix), according to ([75|), is a de- 
creasing function of the hole doping concentration x. The 
underdoped region corresponds to the tightly-bound pre- 
existing-pair regime in this picture; the overdoped region 
is identified with the extended-pair BCS-like regime. The 
superconducting transition temperature Tc as a function 
of x reaches a maximum at an optimal doping concen- 
tration Xo not too far from the concentration xq where 
f^o(a^n) — ^ and the isotope exponent vanishes. Notice 
here that the point Xo = Xo(S^) [where T'c(f2, x) reaches its 
maximum as a function of x at fixed f2] need not exactly 
coincide with the point xq [where the optimal phonon 
frequency rio(a^) equals the actual phonon frequency il]. 

At sufficiently large hole doping the polaron-polaron 
wavefunction overlap and the mutual screening of the 
hole-localizing potential wells Cu may become so strong 
that the holes become unbound, that is, the polarons be- 
come unstable towards forming free carriers. This finite- 
density polaron unbinding can be regarded as analogous 
to the Mott delocalization transition in moderately doped 
semi-conductors. The primary difference is that the Mott 
transition in semi-conductors involves the screening of lo- 
calizing potential wells due to static impurities whereas, 
in the present case, the localizing potential wells are due 
to local lattice distortions which are induced, via the EP 
coupling, by the polaronic holes themselves. In the adia- 



batic potential W{u), this unbinding will manifest itself 
in the (gradual or abrupt) disappearance of local min- 
imum configurations u^^\ Whether, in the thermody- 
namic limit, this occurs as a sharp transition or as a 
continuous cross-over is presently unclear and needs fur- 
ther study [|l2| . The nature of the polaron unbinding and 
the characteristic concentration Xu where the unbinding 
occurs will also be infiuenced by the long-range Coulomb 
interaction Vc and, in more general EP models, by the 
spatial range of the EP interaction, p^ ] 

If the optimal polaronic doping concentrations Xo and 
Xo are close to the polaron unbinding concentration Xu, 
the polaron unbinding will likely dominate the cross- 
over into the extended pair regime: In this scenario 
(xu = Xo,xq), the cross-over from under- to overdoping 
takes the system directly from the tightly-bound polaron 
pair liquid into a BCS-like superconductor of extended 
pairs of non-polaronic carriers. The effective mass of the 
non-polaronic carriers in the overdoped regime is much 
less enhanced by the electron-phonon coupling and, more 
importantly, the mass enhancement is independent of the 
isotopic mass of the ions. The latter is suggested by 
the conventional weak-coupling electron-phonon theory 
where the mass enhancement factor is given by (1 -I- Az) 
and the Eliashberg parameter is independent of the 
isotope mass. If the pairing attraction is of predomi- 
nantly electronic (i.e., non-phonon) origin, one will then 
obtain a very small isotope exponent |a| ^ 1 throughout 
the overdoped regime x > x^. 

Thus, the overall magnitude of a can serve as a distin- 
guishing feature between extended pairs of polaronic and 
non-polaronic carriers in the overdoped regime. In the 
former scenario, already described above, a{x) changes 
sign near optimal doping, but |a| well inside the over- 
doped regime can become as large as in the underdoped 
regime, reflecting the fact that the underlying pair con- 
stituents are still single-hole polarons. By contrast, in 
the latter (unbound carrier) scenario, |q;| becomes small 
in the overdoped regime, without necessarily incurring 
a sign change in a, reflecting the non-polaronic nature 
of the pair constituents. Further experimental studies of 
the isotope exponent in the overdoped cuprate systems 
would be desirable. 

The foregoing features of the underdoped polaron liq- 
uid model and its cross-over into the overdoped regime 
exhibit strong similarities with the observed pairing sym- 
metry and doping dependences of Tc, isotope exponent 
and pseudo-gap in the cuprates. However, in its present 
form, the model also suffers from several potential draw- 
backs which arise from the smaZ/-polaron character of 
the self-localized hole. Small-polaron formation neces- 
sarily implies bandwidths Bi and B2 which cannot be 
substantially larger than the phonon energy scale fi, as 
shown in Figs. p[- |l^. As a consequence, small-polaron 
carriers may be easily localized by disorder and/or long- 
range Coulomb interaction effects. Also, by Eq. ([Tq), the 
overall magnitude of the optimal Tco ~ \xB2 ^ ^xfl 
cannot exceed some fraction of J7. With x ^ 0.10 — 0.20 
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and ^ lOOOX ||], this upper limit on is of order 
50 — lOO-ftT in the cuprates and it is reached if E-p just 
barely exceeds E'p"*^ For substantially larger E-p, B2 
and Tc are rapidly (exponentially !) suppressed with Ep. 
It is not clear from the experimental data, whether ob- 
served carrier mobilities, effective masses and TcS in the 
underdoped cuprates actually exhibit such a strong sen- 
sitivity to changes in Ep and/or to disorder or long-range 
Coulomb interactions. 

The foregoing limitations of the small-polaron system 
can ultimately be traced back to the short-range nature 
of the assumed Holstein EP coupling in our model. Scal- 
ing arguments show that, at the level of the 0th order adi- 
abatic approximation in spatial dimensions D > 2, short- 
range EP models are subject to a dichotomy whereby 
single carriers either form small polarons, if Ep exceeds 



a certain threshold E^ 



(crit) 



larons at all, if Ep < Ep"*"^ 



> 0, or they do not form po- 

||2^ ]. By contrast, in systems 
with additional longer-range EP couplings, such as the 
Frohlich model |^,^, as well as in ID short-range EP 
it is possible to form large polarons at 
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arbitrarily weak Ep, i.e., with iJp'^"*'' = 0. It has been 
argued [ p9| that large-polaron and large-bipolaron mod- 
els can remedy some of the above described deficiencies 
of the small-polaron picture, while retaining most of the 
desirable physical features described above. Thus, in a 
large-polaron model, there is still preformed pair forma- 
tion above the superconducting Tc in the underdoped 
regime; and the carrier bandwidth is still strongly re- 
duced and dependent on isotope mass. Also, the possi- 
bility of a cross-over to a BCS-type free carrier supercon- 
ductor, as a function of doping, via polaron unbinding, is 
retained in a large-polaron theory. However the depen- 
dence of the bandwidth and Tc on EP coupling strengths 
Ep and on phonon frequencies O is only algebraic, rather 
than exponential, and the overall magnitude of the large- 
polaron bandwidths can become substantially larger than 
in a small-polaron model, thus allowing for larger T^'s. 
The proposed large-polaron theories studied so far |3|] 
have been based on phenomenological continuum models 
which of course cannot reproduce lattice- related features, 
such as the location of band minima and pairing sym- 
metries discussed above for our 2D lattice model. It will 
therefore be of interest to extend our present work to lat- 
tice models with longer-range EP couplings. Such future 
studies should explore the possibility of large-polaron 
formation and the bandstructure and pair wavefunction 
symmetry of large-polaron pairs. 

Another critical problem in the above described po- 
laron models is the inclusion of long-range Coulomb ef- 
fects. Rough estimates based on a point charge model 
and measured long-wavelength dielectric constants ^] 
suggest that Vc/t in the cuprates could be as large as 
1 — 2, if only the electronic contribution to the dielec- 
tric screening, that is, only eoo, is taken into account. If 
additional screening from phonons, i.e., eq, is included, 
the estimated Vc/t is reduced to 0.15 — 0.3. The for- 



mer, Vc/t ^ 1 — 2, would be sufficient to completely 
suppress the polaron pairing attraction in a system con- 
taining only two isolated holes, that is, in the limit of 
vanishing polaron density. The latter, Vc/t ~ 0.15 — 0.3, 
may be overcome by the AF-mediated 1st neighbor at- 
traction, but the net attraction strength would still be 
substantially reduced by Vc- [^,^ The suppression of 
extended pairing states, such as (1^2 _y2 pairing, by the 
long-range part of the Coulomb interaction is a common 
problem in all extended pairing models which are cur- 
rently under investigation. Recent studies of the metal- 
lic (in addition to insulating dielectric !) screening of 
the extended Coulomb potential at finite doping density 
p2| , p^ have suggested that the screened Coulomb poten- 
tial becomes substantially reduced, or even attractive, at 
doping concentrations x ~ 0.1 — 0.2. However the forego- 
ing studies are based on weak coupling or diagrammatic 
approaches which do not include polaronic strong-EP ef- 
fects. It therefore remains to be seen whether metallic 
screening in a finite-density polaron liquid will be suffi- 
cient to "rescue" the AF-driven pairing attraction from 
the repulsive long-range Coulomb forces. 

It is also worth re-emphasizing the strong phonon 
contribution to the dielectric screening in the cuprates, as 
evidenced by the large measured dielectric constant ra- 
tio eo/eoo ^ 6 ll^,^. This phonon contribution, which 
acts to reduce long-range Coulomb forces, can be equiv- 
alently regarded as a long-range attraction, mediated by 
long-range (dipolar) EP interactions. This long-range 
EP interaction is in fact a primary agent causing (bi- 
)polaron formation in the above-cited |3^ phenomeno- 
logical large-polaron models. It is therefore quite con- 
ceivable that, in a realistic model of the cuprates, both 
AF- and EP-mediated attractions contribute to the over- 
all pairing potential and that the EP contribution may 
even be the predominant one. 



X. SUMMARY 

In conclusion, we have developed a treatment of po- 
laron tunneling dynamics on the basis of a path integral 
formulation of the adiabatic approximation. The adi- 
abatic treatment of polaron tunneling has been tested 
by comparison to exact numerical results for a two-site 
Holstein system. The break-down of the adiabatic ap- 
proach in the anti-adiabatic regime has been discussed 
and the resulting limitations of applicability for long- 
range polaron tunneling processes in lattice models have 
been identified. Using a combination of path integral, 
many-body tight-binding and exact diagonalization tech- 
niques, we have then explored the Berry phases and ef- 
fective matrix elements for single- and two-polaron tun- 
neling, the two-polaron quasiparticle statistics, effective 
two-polaron interactions, and polaron pairing states in 
the 2D Holstein-t J and Holstein-Hubbard models near 
half filling. The effect of 2nd neighbor electron hy- 
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bridization and long-range Coulomb repulsion has also 
been studied. Due to the AF spin correlations, single- 
polaron hopping is dominated by mtra-sublattice 2nd- 
and 3rd-neighbor processes. These processes are strongly 
enhanced by close proximity of a second polaron. The 
Berry phases imply either d^2_y2- or Px(^yywave pair 
symmetries and effective spin-l/2-fermion quasiparticle 
statistics of dopant-induced polaron carriers. For the 
Holstcin-t J and Holstcin-Hubbard models on the 8-, 10-, 
and 16-site clusters, the dx2-y2-wa.ve state is stable for 
two polarons. The second neighbor hopping Ht^ favors 
the d^2_y2-wa.ve pair for t' > 0, while the long-range 
Coulomb repulsion Hic favor the Px(y)-wa.ve pair. 

The strong on-site Hubbard-C/ Coulomb repulsion 
plays a crucial role in the formation of these pairing 
states. By keeping the electrons spatially separated and 
preventing on-site bi-polaron formation the Hubbard-?/ 
interaction acts, effectively, to greatly enhance the po- 
laron tunneling bandwidths and, hence, their mobility, 
in the nearly ^-filled regime. 

For a hypothetical superconducting polaron pair con- 
densate, our results imply qualitative doping depen- 
dences of the isotope effect, Tc and pseudo-gap which 
are similar to those observed in the cupratcs. Potential 
limitations of the present polaron model, arising from 
the short-range nature of the assumed EP coupling, have 
been pointed out. Further studies to include longer-range 
EP couplings, in combination with extended Coulomb in- 
teractions, have been outlined. 
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FIG. 1. 0th and 1st order contributions Wo and Wi to 
the effective adiabatic lattice potential in the 2-site Holstein 
model with 1 electron. In (a), the electronic groundstate en- 
ergy Wo/Ep and the 1st excited state energy Wq^^/Ep of 
the electronic Hamiltonian Ho{u) are shown as functions of 
u-/up, at M+ = for t/Ep = 0, 0.1, 0.3, 0.5 and 0.7. In (b), 
Wi X Ep/n"^ vs. u-/up is shown at u+ = for t/Ep = 0.1, 
0.3, 0.5 and 0.7. 

FIG. 2. Exact and adiabatic results for the tunneling 
energy splitting, 2tp, between the groundstate and the 1st 
excited electron-phonon eigenstate, plotted as a function of 
Ep/Q, in the 2-site Holstein model with 1 electron, for t — 1 
and several values of Ep, as indicated. 



FIG. 3. (a) One-polaron and (b) two-polaron closed tun- 
neling paths and their Berry phase factors. Black circles in- 
dicate the polaron locations for the initial u-configuration of 
the path. The numbers on the two-polaron exchange paths in 
(b) indicate the order of the single-polaron tunneling steps. 

FIG. 4. Important single-polaron tunneling processes with 
matrix elements tp ' to the v-th neighbor sites for (a) P=l and 
(b) P=2 polaron states on the 2D square lattice. 

FIG. 5. (a) Berry phase contributions from single-polaron 
Ist-neighbor processes in the vicinity of a second, static po- 
laron (black circle). Full and dashed bonds indicate (-1-1) 
and ( — 1) Berry phase contributions, respectively, (b) Inter- 
nal parity of 2nd- and 3rd-neighbor polaron pairs is odd under 
reflection along the dashed line. 

FIG. 6. (a) One-polaron and (b) two-polaron open tun- 
neling paths and their Berry phase factors. Black circles in- 
dicate the polaron locations for the initial u-configuration of 
the path. T{a,b) denotes the translation by vector (a, 6), R{<j)) 
the rotation by angle (f). 

FIG. 7. Total momentum (px,Py) of the (a) P=l and 
(b) P=2 polaron ground state of Hp in a tp-'-vs.-tp-' phase 
diagram. Also shown in (b) are the internal symmetries of 
the respective two-polaron ground states. 

FIG. 8. (a) Logarithm of the effective polaron hop- 
ping amplitudes tp' for P=l and 2 holes and i/=2nA- and 
3rd- neighbor processes vs. inverse phonon energy 1/Q. (with 
t^p' ~ t^p' due to accidental cluster symmetries), (b) Effec- 
tive polaron nearest-neighbor attraction Vp and two-polaron 
binding energy A vs. 1/Q.. All results are for t = 1, J = 0.5t, 
with Ep = C'^/K = 2.5t and 4.0t, on an iV=8 lattice with 
periodic boundary conditions. 

FIG. 9. (a) Logarithm of the effective polaron hop- 
ping amplitudes tp' for P=l and 2 holes and i'=2nd- and 
3rd- neighbor processes vs. inverse phonon energy 1/Q. (with 
ip' = t^p' due to accidental cluster symmetries), (b) Effec- 
tive polaron nearest-neighbor attraction Vp and two-polaron 
binding energy A vs. l/fi. All results are for t = 1, U = 9>t, 
with Ep = 2.5t and 4.0t, on an N=S: lattice with periodic 
boundary conditions. 

FIG. 10. Logarithm of the effective polaron hopping 
amplitudes fj"' for i^=2nd- and 3rd-neighbor processes vs. 
inverse phonon energy with and without inclusion of 

next-nearest-neighbor hopping t'=~\-Q.2t, for t = 1, J — 0.5t, 
(a) Ep = 2.5t, and (b) Ep = 4.0t, on an N=8 lattice with 
periodic boundary conditions. 
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FIG. 11. Logarithm of the effective polaxon hopping 
ampHtudes tj'^' for j/=2nd- and 3rd-neighbor processes vs. 

inverse phonon energy l/f2, with and without inclusion of 
next-nearcst-neighbor hopping t'=—0.2t, for t = 1, J = 0.5t, 
(a) Ep = 2.5t, and (b) Ep = A.Ot, on an N=8 lattice with 
periodic boundary conditions. 

FIG. 12. Logarithm of the effective polaron hopping 
amplitudes t^^ for i^=2nd- and 3rd-neighbor processes vs. 
inverse phonon energy with and without inclusion of 

long-range repulsion Vc = l-Ot, for t = 1, J = 0.5t, (a) 
-Ep — 2.5t, and (b) Ep — 4.0t, on an N=8 lattice with peri- 
odic boundary conditions. 
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